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VECTOR  POLAR  METHOD  FOR  THE  EVALUATION 
OF  WAVE  INTERACTION  PROCESSES 

Abstract 

Vector  polar  method  a  graphical  technique  for  the  solution  of  wave  interaction  problems 
that  are  concerned  with  the  overall  results  rather  than  with  the  progress  of  the  process. 
The  method  is  based  on  the  use  of  wave  polars,  that  is  diagrams  representing  the  l(x:i  of 
states  attained  by  a  given  wave  process  in  a  coordinate  system  of  which  one  component 
is  pressure  or  velocity  of  sound  and  the  other  is  particle  velocity.  Furthermore,  if  the  ratio 
of  pressures  or  sound  speeds  is  represented  on  a  logarithmic  scale,  usually  as  the  ordinate, 
while  the  abscissa  represents  the  change  in  particle  velocity  across  the  wave  in  a  linear  scale, 
then  such  a  generalired  hodograph  plane  has  the  property  of  a  vector  field,  that  is,  the  change 
of  state  brought  about  by  the  wave  action  is  evpressed  in  terms  of  a  vector  between  the 
origin  and  any  point  of  the  polar.  By  simple  vector  addition  one  can  then  evaluate  rapidly 
some  quite  complex  wave  interaction  processes. 

The  utility  of  the  method  is  illustrated  here  by  a  variety  of  shock-tube  problems  involv  ing 
interactions  between  sluxtks,  deflagrations,  rarefaction  waves,  contact  discontinuities  and 
area  change.  Although  the  analysis  of  most  of  such  problems  is  known,  the  application* 
of  the  vector  polar  method  led  to  the  discovery  of  some  interesting  regimes  of  solutions  that 
would  be  difficult  to  discern  by  the  more  laborious  conventional  techniques. 
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I.  Introduction 

One  of  the  most  useful  ways  of  treating  problems  in  non-steady  gas  dynamics  is  to 
consider  them  as  results  of  wave  interaction  processes.  A  wave  is  a  regime  that  propagates 
through  the  medium  producing  a  change  in  state,  subject  to  restrictions  imposed  by  the 
principles  of  conservation  of  mass,  species,  momentum  and  energy.  If.  in  particular,  one 
wishes  to  know  only  what  stationary  states  are  created  by  a  given  wave  interaction  pro- 
ccjS  without  inquiring  into  the  details  of  its  progress,  one  should  find  the  method  intro¬ 
duced  here  particularly  useful. 

It  is  sufficient  for  this  purpose  to  consider  at  first  two  elementary  kinds  of  wave 
processes:  one  of  the  steady,  or  ”Hugoniot-typc”,  and  the  other  non-steady,  or  ..Ric- 
mann’ian’*.  The  first  is  distinguished  by  the  fact  that  wave  characteristics  in  the  time- 
space  domain  (wave  world-lines)  are  parallel,  and  includes  the  classical  case  of  the 
gasdynamic  discontinuity;  the  second  is  represented  in  its  most  elementary  form  by  the 
simple,  or  uniform  wave  whose  characteristics  are  straight  lines. 
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Fundamental  properties  of  these  waves  are  derived  here  from  basic  principles,  and  the 
change  of  state  brought  about  by  their  action  is  described  in  terms  of  polar  diagrams 
which  are  then  used  for  the  vector  polar  method  —  a  graphical  technique  for  the  evaluation 
of  wave  interaction  processes.  In  order  to  treat  wave  interactions  with  area  changes,  the 
process  of  steady-state  flow  had  to  be  also  represented  by  means  of  suitable  polar  diagrams. 

Foundations  for  the  finite,  wave  interaction  analysis  were  laid  down  by  Riemann  (1859), 
Hugoniot  (1887—9)  and  Hadamard  (1903).  After  some  early  attempts  at  a  step-wise 
solution  of  some  non-steady  gasdynamic  problems,  such  as  those  of  Aschenbrenner  (1937) 
and  Pfriem  (1941),  the  analysis  has  been  reduced  to  an  engineering  reality  by  such  con¬ 
tributions  as  those  of  Schultz-Grunow  (1944),  Sauer  (1942),  de  Haller  (1945),  Jenny 
(1949)  and  many  others  who  in  the  1940’s  have  popularized  for  this  purpose  the  use  of 
the  method  of  characteristics.  Since  then  such  considerations  became  the  nucleus  of  the 
rapidly  expanding  field  of  shock  tube  technology  with  publications  too  numerous  to 
mention  here.  Most  sicnificant  to  the  subject  matter  of  this  paper  are  the  comprehensive 
texts  of  Courant  and  Friedrisch  (1948),  Rudinger  (1955),  Foa  (I960),  Glass  and  Hall 
(1958),  Oswatitsch  (1956),  Von  Mises  (1958),  Stanyukovich  (1955),  and  Wright  (1961), 
each  containing  quite  an  exhaustive  bibliography.  The  importance  to  engineering  practice 
of  the  non-steady  flow  phenomena  that  can  be  treated  by  our  methods  is  demonstrated 
by  the  papers  of  Bannister  and  Mucklow  (1948),  Wallace  and  Mitchell  (1952)  and  Mucklow 
and  Wilson  (1955). 

The  application  of  the  method  of  characteristics  to  the  solution  of  wave  interaction 
processes  is  restneted  in  scope  to  the  continuous  domain  of  the  phenomenon.  The  method 
can  be  therefore  used  only  to  trace  the  details  of  the  process.  In  contrast  to  such  an  ana- 
^  lysis  ”in  the  small”,  it  is  often  of  interest  to  know  the  results  of  the  interaction  without 
tracing  the  details  of  its  progress,  in  other  words  to  evaluate  the  stationary  states  and  the 
.wave  system  that  has  been  attained  by  a  given  wave  interaction  after  all  the  transient 
phenomena  have  died  down.  Moreover  the  answer  t;^  such  a  question  yields  information 
on  the  final  bound  to  the  solution  of  the  inherently  initial  value  problem  of  the  hyperbolic 
domain  treated  by  the  method  of  characteristics. 

It  is  for  the  purpose  of  getting  such  an  answer  directly  that  the  vector  polar  methods 
is  most  suitable.  In  contrast  to  the  method  of  characteristics,  it  provides  only  an  analysis 
”in  the  large”,  leaving  the  zone  of  interaction  completely  outside  the  scope  of  inquiry. 
In  many  practical  instances  however  this  is  quite  sufficient,  especially  if  the  duration  of 
the  wave  interaction  event  is  short  in  comparison  to  the  rest  of  the  wave  propagation 
processes  under  consideration.  In  our  experience  this  was  indeed  true  in  the  case  of  wave 
interaction  phenomena  that  occur  in  the  course  of  the  development  of  detonation  in 
gaseous  mixtures  [Oppenheim  and  Stern  (1958),  Oppenheim,  Stern  and  Urtiew  (I960)]. 

As  one  would  suspect,  there  exists  an  equivalent  method  for  the  evaluation  of  steady- 
flow,  two-dimensional  wave  interaction  processes  (see  e..g  Fallis  et  a!  (1953));  however 
its  usefulness  is  handicapped  by  the  awkwardness  of  polar  dependence  on  initial  Mach 
number  and  the  inherent  ambiguity  associated  with  the  possibility  of  a  subsonic  or  a  super¬ 
sonic  solution  that  can  be  obtained  from  the  oblique  shock  polar  for  a  given  deflection 
angle. 

With  respect  to  non-steady  flow,  it  is  of  interest  to  note  that  a  method  quite  similar 
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to  ours,  but  more  restricted  in  scope,  has  been  published  independently  in  Poland  by 
Kijkowski  (1959). 

For  the  sake  of  an  orderly  presentation  the  polars  are  here  introduced  first,  and  then 
their  application  to  various  wave  interaction  processes  is  illustrated  by  several  examples. 
It  should  be  noted  however  that  the  paper  does  not  have  to  be  read  in  the  same  sequence 
us  it  has  been  written.  In  fact,  the  reader  may  find  it  more  amusing  to  acquaint  himself 
f  rst  with  solutions  that  may  strike  his  interest,  before  learning  about  the  derivation  of 
the  polars  and  their  properties. 


Notation 


a  —  local  velocity  of  sound, 

A  —  velocity  of  sound  ratio 

—  specific  heat  at  constant  volume, 
e  —  internal  energy, 
h  —  enthalpy, 

L  —  reference  length, 

Af  —  Mach  number, 
p  —  pressure, 

P  —  pressure  ratio  (p,IPx), 
q  —heat  of  reaction, 

Q  —  non-dimensional  heat  of  reaction  (qlp^v^, 

R  —area  ratio, 
t  —  time, 

T  -  non-dimensional  time  coordinate  OoijL, 

T  —  absolute  temperature, 

M  —  particle  velocity, 

U  —  non-dimensional  particle  velocity  [  iu/Ox  =  (m,,  —  m,)  a,], 
r  —  non-dimensional  particle  velocity  (Jm/| 

V  —  specific  volume, 

)  ’  —  specific  volume  ratio, 

V  —  wave  velocity  relative  to  medium  into  which  it  propagates, 
y  —  non-dimensional  wave  velocity  (v/a,), 

F  — non-dimensional  wave  velocity  (v/( 

If  —  absolute  wave  velocity  («  f  a), 

\y  —  non-dimensional  absolute  wave  velocity  {U+A), 

X  —  distance, 

A'  —  non-dimensior  1  sp.ace  coordinate  (.v/L), 

7  —  specific  heat  ratio, 

0  -  density. 


Subscripts 

a  —  refers  to  the  initial  velocity  of  sound, 
c  —  denotes  convergence. 
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—  denotes  divergence, 

—  denotes  incident, 

—  refers  to  reflected  shock, 

—  refers  to  a  shock  front  (i.c.  M,  its  Mach  number  with  respect  to  state  into 
which  it  propagates,  P,  pressure  ratio  across  the  wave), 

—  refers  to  the  initial  state  at  rest, 

—  refers  to  state  immediately  in  front  of  wave  or  initial  state  for  a  process, 

—  refers  to  state  immediately  behind  the  wave  or  final  state  for  a  process, 

.  —  denote  states  indicated  on  diagrams, 

—  denotes  critical  value, 

—  denotes  throat  secti  on. 


2.  Wave  PoUrs 

2.1.  Steady  Wave 


A  Steady  wave  is  one  that  propagates  through  space  at  a  constant  velocity  without 
any  distortion  in  its  structure.  Physically  it  is  a  plane  wave  whose  world-lines  (traces  in 
space-time)  are  loci  of  constant  states  forming  a  family  of  straight  parallel  lines,  as  shown 
in  Fig.  1.  Notable  examples  of  such  waves  are  shocks,  deflagrations  and  detonations. 


Fig.  1.  Steady  Wave  in  the  Space-Time  Domain 

Basic  wave  relations,  namely  the  expressions  for  the  change  in  particle  velocity  that 
it  produces,  and  for  the  velocity  of  its  propagation,  arc  obtained  directly  from  the  conti¬ 
nuity  and  momentum  equations  which,  with  reference  to  Fig.  1,  can  be  written  respectively 
as  follows: 

-Jv 

V 


(1) 
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and 


(2)  ~Ju=Jp 

where  =  m,  —  Uj,,  Jv  =  v,  —  Vj,  and  Ap  =  p^—  />,. 

The  two  relations  are  then  obtained  directly  by  multiplication  and  division,  yielding: 

(3)  (duY  =  —  Ap‘  Av 
and 


(4) 


-V* 


Av 


The  appropriate  Hugoniot  equation  for  the  wave  process  gives  a  functional  depen¬ 
dence  Ap  =  f{Av)  with  the  use  of  which  (3)  and  (4)  become  functions  of  a  single  variable, 
yielding  the  relation  for  the  polar:  Au  =  u{p)  and  for  the  propagation  speed  v,  =  v(p). 

It  should  be  noted  that  expressions  (3)  and  (4)  are  walid  irrespectively  of  the  thermo¬ 
dynamic  properties  of  the  substance,  and,  within  the  dynamic  meaning  of  p  as  the  local 
normal  component  of  the  stress  tensor  and  of  v  as  the  reciprocal  of  local  density,  they 
apply  without  invoking  any  principle  of  equilibrium.  The  most  unambiguous  non-dimen- 
sionalizing  parameter  for  such  a  general  case  is  the  product  p^v^  referred  to  the  state  of  the 
medium  into  which  the  wave  propagates.  Introducing  then  U  =  Aul^ pjo^y  V  = 

P  =  Pi\P*  and  V  =  Vy/Vg,  equations  (3)  and  (4)  become: 

(3a)  (/*  =  (P-I)(1-T;) 

(4a)  K*=  (P- I)/(I -Q>) 

For  a  perfect  gas  with  constant  specific  heats  the  Hugoniot  equation  can  be  written  as 
follows: 


(5)  — L  pv--t.--p.v,-q=  (<;  +  »,) 

y- I  r.- 1  2 

where  q  is  the  heat  released  per  unit  mass.  With  Q  =  qlp^v^,  it  becomes: 


(5a) 

which  in  terms  of 

and 


y-  1 


y.-i 


reduces  to: 
(5b) 


(/>  +  ^)(i;-^)  =  (5  +  ^)(i-^) 


demonstrating  that  the  Hugoniot  curve  is  a  simple  hyperbola  with  asymptotes 'T/’ =  ^ 
and  /»=  —  /?. 
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Equations  (3a)  and  (4a)  yield  now,  with  the  use  of  equation  (5b),  the  following  expres¬ 
sions  respectively  for  the  wave  polar  and  for  the  wave  propagation  velocity: 

(6) 

iP  +  ^) 

(7)  + 

which  summarize  all  the  salient  properties  of  these  curves.  For  example,  the  condition 
of  existence:  P  >  B  P  >  I  which  is  evident  from  the  form  of  (6)  and  (7)  corresponds 
to  the  well  known  exclusion  of  the  ’’forbidden  zone”  that  is  clearly  delineated  by  (5b) 
as  one  implying  that  P  >  B  can  be  attained  only  if  <  1. 

It  is  customary  for  a  perfect  gas  to  express  velocities  in  non-dimensional  form  by  re¬ 
ference  to  the  velocity  of  sound.  Equations  (6)  and  (7)  give  immediately  the  required 
expressions  for  the  change  in  particle  velocity: 


and  for  the  wave  Mach  number: 


I  A/£--=— r- 

0*  y  Yx 

The  ratio  of  the  local  velocities  of  sound  is: 


or,  with  (5b): 


Yx  I  (P  +  P)  J  Yx 


Yx  Y+  1 
Y-l 


B^P+\ 


If  Yx  =  Yp  then  =  1  +  (y  —  1)^  =  1  +  j — ^  j.  For  a  shock  wave  ^  =  0,  hence 
B  =  1  and  equation  (6)  [with  equation  (8)]  becomes: 

(11)  (/*  =  __  2  _ 

y(y  +1)  I) 

Y  +  1 

while  equation  (7)  [with  (9)]  reduces  to  the  well  known  relation  between  the  Mach  number 
and  pressure  ratio  for  a  normal  shock  wave 


P-  \ 

y+  1 
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Fig.  2.  Shock  Polars  for  Perfect  Gases  in  the  Pressure  Hodograph  Plane  and  Shock  Mach  Numbers  as 

a  Function  of  Pressure  Ratio 


At  the  same  time  equation  (10a)  yields  a  simple  relation: 


(13) 


/(*  = 


Steady  wave  polars  are  introduced  in  Figs.  2  to  7.  With  a  logarithmic  scale  for  the  pre¬ 
ssure  ratio  and  for  the  speed  of  sound  ratio,  while  a  linear  scale  is  maintained  for  the  velo¬ 
city  change  across  the  wave,  the  polar  diagrams  are  rendered  a  vectorial  character.  They 
form,  in  fact,  generalized  velocity  hodographs,  each  polar  representing  a  locus  of  end 
points  of  vectors  starting  from  the  origin  whose  one  component  is  the  change  of  particle 
velocity  across  a  given  wave  V  =  /Suja,  and  the  other  is  the  logarithm  of  pressure  ratio 
P  =  PylPx  (or  the  sound  speed  ratio  A  =  fl,/o,). 

Fig.  2  is  a  plot  of  equations  (1 1)  .':nd  (12)  for  various  values  of  specific  heat  ratios,  y, 
from  1  to  5/3.  Fig.  3  is  a  corresponding  plot  in  the  domain  of  the  local  velocity  of  sound. 


I 


448 


A.  K.  OPPENHEIM,  P.  A.  URTIEW  AND  A.  J,  LADERMAN 


HI 


— d 


Fig.  3.  Shock  Polars  for  Perfect  Gases  in  the  Velocity  of  Fig.  4.  Shock  Polars  and  Mach  Number  for  a  Persfect  Gas 

Sound  Hodograph  Plane  with  y  =  1.4  in  the  Pressure  Hodograph  Plane 
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Fig.  5.  Shock  Polars  for  a  Perfect  Gas  with  y  =  1.4  in  the  Velocity  of  Sound  Hodograph  Plane 


Ua  '/o,  -  — 


Fig.  6.  Deflagration  Polars  and  Mach  Number  for  the  Case  od  =  1.4,  =  1.2  :  dQ  =  42.28  (repre¬ 

senting  an  H,  + 1  2  O,  mixture)  in  the  Pressure  Hodograph  Plane 

The  line  of  A/  =  I,  which  in  Fig.  3  is  simply  a  logarithmic  curve,  marks  the  limiting 
strength  of  shocks  below  which  the  flow,  with  respect  to  observer  who  is  stationary  in  the 
regime  into  which  the  wave  propagates,  is  subsonic,  and  above  which  it  is  supersonic. 

Figs.  4  and  5  are  the  working  diagrams  for  the  solution  of  problems  involving  wave 
interaction  processes.  In  this  case  the  velocity  change  has  to  be  related  to  some  constant 
reference  value  of  the  velocity  of  sound  which  is  usually  taken  as  that  of  the  undisturbed 

Au 

medium  initially  at  rest.  Consequently  the  abscissa  in  the  working  diagrams  is  t/,  =  = 
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Fig.  7.  Deflagration  Polars  for  the  Case  of  =  1 .4,  y,  =  1 .2  and  Q  «=  42.28  in  the  Velocity  of  Sound 

Hodograph  Plane 

=  -  -  -  =  UAg  and  each  polar  becomes  a  function  of  the  velocity  of  sound  in  the  medium 

Ox  Oq 

into  which  the  wave  propagates.  Figs.  4  and  5  represent  then  such  a  family  of  polars  de¬ 
rived  quite  simply  from  those  of  Figs.  2  and  3  for  y  =  1.4  and  /I,  of  1,  2,  3,  and  4  with 
interpolation  scales  for  fractional  values  of  Ax- 

Figs.  6  and  7  represent  working  polar  diagrams  for  a  deflagration  wave  propagating 
to  the  right.  The  left  hand  side  of  Fig.  6  is  essentially  the  plot  of  equation  (6)  [with  (8)] 
the  right  hand  side  is  the  plot  of  equation  (7)  [with  (9)]  and  Fig.  7  is  transformed  from 
Fig.  6  by  the  use  of  equation  (10a).  Under  the  assumption  of  perfect  gas  behavior  of 
reactants  and  products,  the  deflagration  is  completely  defined  in  terms  of  y,  =  1.4,  = 

=  1.2,  and  Q  —  42.3  that  correspond  here  to  the  properties  of  a  stoichiometric  hydrogen- 
oxygen  mixture,  the  value  of  heat  release  having  been  adjusted  so  that  the  polar  curves 
would  reproduce  the  experimental  value  of  the  Chapman-Jouguet  detonation  velocity. 


2.2.  Simple  Wiee 

A  simple  wave  is  characterized  by  the  fact  that  its  world-lines  are  loci  of  constant 
states  forming  a  family  of  straight  (but  not  parallel)  lines  as  shown  in  Fig.  4  which  illu¬ 
strates  two  intersecting  simple  waves.  The  regime  of  interaction  is  an  example  of  an  unsteady 
flow  process  that  is  outside  the  scope  of  our  inquiry. 
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The  two  basic  relations  for  the  simple  wave  can  be  obtained  by  differentiation  from 
those  for  a  steady  wave.  Thus  the  process  relation  is  from  (1) 


(14) 

or  from  (2) 


,  dv  do 

du  =  —  a  -=  a—' 


V 


du  = 


vdp 


while  the  expression  for  the  local  wave  velocity  is  from  (4); 
(15) 


^2  _  _  ,.2 


dv  dn 

Again  the  simplest  non-dimensionalizing  parameter  is  the  local  pv  product  so  that 
(14)  and  (15)  become: 

du  ^/(log  p) 
pv  a 


(14a) 

and 

(15a) 


P^' 


vdp 

pdv 


vdp 


r: 


In  the  integral  form  the  process  relation  is,  from  (14); 

« 

(16)  Ju  = 

1 

giving  the  generalized  von  Mises  expression  for  the  Riemann  relation. 
For  an  isentropic  process  the  above  becomes: 

s 

(17)  Ju  = 

t 


while  (1 5a)  yields; 
(18) 


a 


dh 


pv  de  e 


dWe) 

de 


For  a  perfect  gas,  for  which  — =  0,  we  obtain  then  the  familiar  expression 


de 


(19) 
or 

so  that  (17)  reduces  ti 

(20) 


pv  e 


a'^  =  ^  h  =  (}'  —  \)h 


Au  = 


lada 


(y— l)a  y-1 


Aa 


y 

! 
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The  above  is  the  well  known  Riemann  invariance  w  hich  in  the  present  case  represents 
the  process  relation  for  the  wave. 

The  expression  for  the  local  velocity  of  wave  propagation  in  terms  of  pressure  is  obtained 
directly  from  the  iscntropic  relation  for  the  ratio  of  local  velocities  of  sound 

a 

(21)  A=  ^  = 

Ox 

The  process  relation  (20)  becomes  then; 

\ii  2  "> 

(22)  I)  ‘  ^p  t,  _  1) 

o,  >'  -  1  j*  I 


A  ♦ 

Fig.  8.  Simple  Waves  in  the  Space-Time  Domain 


Fig.  9.  Rarefaction  Polars  for  Perfect  Gases  in  the  Pressure  Hodograph  Plane  and  Local  Velocities  of  Sound 

as  a  Function  of  Pressure  Ratio 
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Fig.  9  represents  the  plot  of  equation  (21)  on  the  left  side,  and  that  of  equation  (22) 
on  the  right  side  for  a  simple  trefaction  wave  propagating  to  the  right.  Fig.  10  is  the 
corresponding  polar  diagram  in  the  domain  of  the  velocity  of  sound,  representing  essen¬ 
tially  a  cross-plot  of  the  curves  of  Fig.  9.  Similarly  as  in  Figs.  2  and  3,  the  family  of  polars 
are  given  here  in  terms  of  the  specific  heat  ratio  as  the  parameter,  covering  the  whole 
range  of  y  from  I  to  5/3. 

f-igs.  11  and  12  are  working  diagrams  obtained  from  Figs.  9  and  10  for  y  =  1.4  in 
exactly  the  same  way  as  Figs.  4  and  5  were  derived  from  Figs.  2  and  3. 


2.3.  Steady  Flow  Polars 

In  order  to  apply  our  graphical  method  to  problems  involving  wave  propagation 
through  ducts  with  changes  in  cross-section  area,  the  steady  flow  has  to  be  specified  in 
terms  of  polar  diagrams  as  well.  Such  diagrams  are  represented  in  Figs.  13  to  16  for  an 
iscntropic  flow  of  a  perfect  gas  witli  y  =  1 .4.  The  meaning  of  the  polars  is  the  same  as 


Fig.  13.  Steady  Flow  Polars  for  Divergences  in  the  Pressure  Fludograph  Plane  (y  —  1.4) 
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before,  that  is  they  represent  generalized  velocity  hodographs.  The  scales  are  identical, 
in  order  to  permit  vectorial  superposition.  However,  unlike  the  non-steady  polars,  the 
end  point  of  the  steady  flow  velocity  vector  is  not  only  subject  to  compatibility  conditions 
between  waves  of  the  system  resulting  from  a  given  interaction  process,  but  it  is  also  fi.xed 
by  the  value  of  the  Mach  number  at  the  inlet.  This,  on  one  hand,  gives  rise  to  the  genera¬ 
tion  of  new  waves,  and  on  the  other,  necessitates  the  use  of  auxiliary  graphs,  which  have 
been  included  on  the  right  sides  of  Figs.  13  to  16. 


Fig.  16.  Steady  Flow  Polars  for  Convergences  in  the  Velocity  of  Sound  Hodograph  Plane  (•/  1.4) 


The  graphical  procedure  is  illustrated  by  examples  which,  in  Figs.  13  and  14,  demon¬ 
strate  how  to  obtain  the  end  point  of  the  steady  flow  vector  for  the  area  ratio  R  2.5, 
and  the  inlet  Mach  number,  3/,  =  1.5,  while  in  Figs.  15  and  16  for  the  area  ratio  R  0.8 
and  the  inlet  Mach  number  3/j  =  0.5.  As  it  is  evident  there,  the  area  ratio  is  used  first 
to  locate  the  appropriate  polar,  and  the  Mach  number  at  inlet  specifies  then  the  end  point. 

In  terms  of  the  pressure  ratio  he  equation  for  the  steady  flow  polar  is: 


(23)  U  = 

and  fo.  the  auxiliary  curve: 


I  -  (P./rJ 


(24) 
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In  terms  of  the  velocity  of  sound  ratio  the  corresponding  equations  are: 

(25) 
and 


(26) 
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The  specific  features  of  the  steady  flow  polars  will  now  be  described  for  the  subsonic 
and  supersonic  case  respectively. 

All  subsonic  polars  pass  through  the  origin  of  the  coordinate  system  on  the  hodograph 
plane.  On  the  other  side  all  hodographs  are  limited  by  the  condition  of  Mach  number 
of  unity  obtained  simply  from  equations  (23)  to  (26)  for  A/,  —  1  which  is  represented 
in  our  diagrams  by  a  broken  line.  For  the  divergence  the  limiting  polar  is  that  for  R  —  Vj, 
corresponding  to  the  case  of  an  open  end.  The  end  point  of  this  polar  represents  the  state 
of  vacuum  attained  at  the  section  of  infinite  divergence,  M^  0,  for  which  the  pressure 

ratio  and  the  velocity  of  sound  ratio  at  (he  end  of  the  tube,  acting  now  as  the  throat 

/ 

of  the  nozzle,  attain  their  critical  values  2  'j  (~2^)  while  the 


maximum  change  in  particle  velocity  is 


lu 


Oi 


—  f.  In  the  case  of  flow  through  a  conver¬ 


gence  the  limiting  polar  is  that  corresponding  to  R  =  0  and  it  represents  the  case  of 
inflow  from  infinite  reservoir  while  the  coordinates  of  its  end  point  are  the  critical  values: 


P*  - 


/  2  ' 

y-i  / 

\y  '  1 , 

11 
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and  — 
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In  contrast  to  the  subsonic  case,  the  supersonic  polars  have  their  singular  points  at 
the  axis  of  ordinates  U  =  0,  corresponding  to  the  condition  of  A/,  =  x.  The  position 
of  each  polar  at  this  axis  depends  only  on  the  value  of  the  area  ratio,  namely: 


(27) 
and 

(28) 


Lim  P  =  /?■>’ 


Lim  A  =  R 


y-1 

2 


which  are  obtained  directly  from  equations  (23)  and  (25)  for  t'  =  0  and  A/,  ^00. 

The  choking  coiidition  for  steady  supersonic  flow  (i.e.  A/,  =  I.O  in  the  divergent  case 
and  M,=  1.0  in  the  convergent  case)  is  represented  by  the  broken  lines  in  Figs.  13—16. 
For  the  divergent  case  u  approaches  an  asymptotic  value,  which  in  terms  of  the  specific 
heat  ratio  is  expressed  as 


2* 


f 
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(29) 


For  }'  =  1.4  this  results  in  Lim  V  =  1.45.  Asymptotic  values  for  other  specific  heat  ratios 

R-rxj 

My-CO 

are  listed  in  Table  1.  The  choking  line  in  the  convergent  case  has  no  asymptote  since 
both  U  and  P  go  to  infinity. 


Table  I 

Functions  of  Specific  Heat  Ratio 


^  Quantity 

Equ¬ 

1 

Ratio  of  Specific  Heats,  */ 

ation 

1.0 

l.l 

1.2 

1  '-3 

1 

Lim  U 

R  -.00 

29 

cx) 

3.58 

2.32 

1.77 

1.45 

1 

1.00 

t 

39 

1.619 

1.697 

1.797 

1.918 

2.069 

2.757 

Lim  A/( 

R->oa 

38 

1.0 

1  1.09 

1.125 

1.145 

1 

I.I5I 

1.17 

1 

Lim  R^ll 
i  W|  ->(xi 

46 

0 

0  008 

0.13 

0.375 

1 

0.648 

0.933 

Lim  Riv 

-»ou 

47 

0 

0  516 

0.68 

0.772 

0.838 

0.957 

Lil  1  A/( 

R  -.00 

— 

1.271 

1.37 

1.47 

1.58 

1.719 

2.075 

Lim  R||i 

Mi  -.00 

57 

00 

12.5 

7.7 

2.67 

1.543 

1.072 

P. 

60 

2.62 

7.28 

8.7 

9.43 

10.35  ' 

12.75 

Lim  (F,),v 

M  — 

65 

2  615 

4.7 

5.12 

5.54 

5.96 

7.16 

3.  Examples 

The  vector  polar  methods  is  illustrated  here  by  the  following  elementary  exaples:  the 
shock  tube  problem,  the  interaction  between  a  shock  wave  and  a  deflagration  front,  the 
merging  of  shocks  and  the  merging  of  a  rarefaction  wave  with  a  shock  front. 


3.1.  Thr  Shock  Tobv  Problem 

The  simplest  case  of  wave  interaction  is  that  of  the  start  of  a  simple  shock  tube.  Before 
the  interaction  there  arc  two  regimes,  both  at  rest,  each  at  a  different  pressure,  separated 
by  a  diaphragm.  The  diaphragm  is  suddenly  removed,  causing  a  shock  to  move  into  the 
low  pressure  and  a  rarefaction  front  into  the  high  pressure  regime. 
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In  the  example  of  Fig.  17  the  initial  diaphragm  pressure  ratio  was  ■  -  2.62.  The 
medium  on  both  sides  of  the  diaphragm  v  as  kept  at  the  same  initial  temperature,  and  it 
was  assumed  to  behave  essentially  as  a  perfect  gas  with  y  —  1.4.  With  reference  to  the 
solution  depicted  on  Fig.  17  the  fact  that  there  is  a  shock  propagating  into  the  regime 
at  state  0  is  reflected  on  the  hodograph  plane  by  the  shock  polar  for  ^  \  of  Fig.  6 
starting  from  point  L’  =  f\,P  1  and  A  =  I.  At  the  same  time  the  fact  that  a  rarefaction 


Fig.  17.  Simple  Shock  Tube  Problem 


is  moving  into  state  4  is  represented  on  the  hodograph  plane  by  the  rarefaction  polar 
for  /I,  1  of  Fig.  11  starting  from  point  6  =0,  P  =  2.62  and  A  =  \.  The  intersection 

of  the  two  polars  on  the  plane  deterrpines  states  (1,  la)  while  the  corresponding  points 
in  the  A  -  V  plane  specify  the  conditions  on  the  two  sides  of  the  contact  discontinuity 
as  shown  in  Fig.  17  and  in  Table  2.  Included  there  moreover  is  the  reilection  of  the  shock 
wave  from  the  closed  end,  establishing  state  2,  and  the  reflection  of  the  rarefaction  fan 
from  the  other  side  producing  state  3.  These  two  states  are  determined  directly  by  the 
intersection  with  the  V  ^  0  axis  of  a  shock  polar  starting  from  state  1  for  Ai  =  1.067, 
and  of  a  rarefaction  wa\c  starting  from  state  la  for  A^^  =  .933. 

The  above  problem  is  indeed  so  simple  that  for  a  perfect  gas  with  constant  specific 
heats  one  can  solve  it  alge>^raically  without  any  difficulty.  In  fact,  according  to  equation  (11) 
the  particle  velocity  behind  the  shock  is  given  by: 
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while  behind  the  rarefaction  wave  equation  (22)  gives 


Since  m,  -  i/,,  and  Pi  =  />i,  one  obtains  then  directly: 


(32) 

w  hence : 

(33) 
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Equation  (26)  yields  the  diaphragm  pressure  ratio,  P^,  that  is  required  to  produce 
a  given  shock  pressure  ratio.  P,,  which  is  related  in  turn  to  the  incident  shock  Mach  number 
by  equation  (12). 

Similarly,  it  can  be  show  n  that  as  a  consequence  of  the  requirement  that  particle  velo¬ 
city  at  the  closed  end  of  the  tube  must  be  zero  behind  the  reflected  shock  wave, 


37-1 


(34) 


{2ft  ^  \)P,-ft  _  7 

Pt  —  Pi  ~  Pi  ~ 

ftp,  f  I 


Pi-1 

7^- r 

r-  1 


yielding  a  straight  forward  expression  for  the  pressure  attained  behind  the  reflected 
wave  in  terms  of  a  given  pressure  ratio  of  the  incident  shock. 

It  should  be  noted  that  the  algebraic  solution  follows  in  essence  the  graphical  treatment. 
The  latter  moreover  has  the  advantage  of  simplicity  and  generality  since,  as  contrasted 
to  the  former,  it  can  be  applied  as  easily  to  the  inverse  as  to  the  direct  problem,  and  it  is 
not  restricted  in  scope  to  any  particular  form  of  the  equation  of  state. 


3.2.  Inlertctioa  Between  a  Shock  and  a  Dellairitioa  Waee 

Consider  a  head-on  collision  between  a  shock  moving  to  the  right  at  A/,  =  1.6  and 
a  deflagration  moving  to  the  left  with  M  =  0.14.  As  shown  in  Table  2  the  shocK  raises 
the  pic.>3i.'re  by  a  factor  of  2.75  while  the  deflagra’ion  is  associated  with  a  pressure  drop 
of  40"o  The  results  of  the  interaction  in  terms  of  stationary  states  that  should  be  attained 
provided  that  the  deflagration  survives  the  shock  impact  arc  showh  in  Fig.  18.  As  it  ap¬ 
pears  there  the  transmitted  deflagration  polar  reaches  the  Chapman-Jouguet  state  (2) 
before  intersecting  the  polar  (10  -20)  of  the  transmited  shock.  The  onl>  way  then  that 
a  compatible  w;  vc  system  would  be  obtained  is  by  postulating  that  the  interaction  is 
associated  with  the  generation  of  a  rarefaction  wave  (2  20)  that  propagates  immediately 

behind  the  Chapman-Jouguet  state. 
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Fif  18.  Interaction  between  a  Shock  and  a  Deflagration  NVave 


The  graphical  technique  sugge  sts  here  directly  the  appearance  of  the  rarefaction  wave 
whose  existence  could  not  have  been  so  easily  deduced  by  analysis. 


3.3.  Skock 

Two  shocks  propagating  in  the  same  direction  must  merge  and,  as  it  has  been  proven 
by  von  Neumann  (1943),  such  an  interaction  results  in  a  transmitted  shock  and  a  reflected 
rarefaction.  As  illustrated  in  Fig.  19,  the  appearance  of  the  rarefaction  wave  is  a  con¬ 
sequence  of  the  fact  that,  for  shocks  propagating  to  the  right,  the  pressure  velocity  ho- 
dograph  has  a  property  that  a  shock  polar  started  from  any  point  of  another  sh-»ck  polar 
must  lie  to  the  left  of  it.  Consequently,  the  only  way  to  attain  a  compatible  wave  system 
is  by  postulating  the  existence  of  a  rarefaction  such  as  2  3  in  Fig.  19.  It  is  of  interest 
to  note  that  the  opposite  holds  true  for  the  sound  speed  velocity  hodograph,  so  that  state  2 
has  a  .3wcr  sound  speed  than  state  3.  Numerical  results  of  our  graphical  solution  are 
again  presented  in  Table  2. 


462 


A.  K  OPPKNHEIM,  P,  A.  URTIEW  AND  A.  J.  LADERMAN 


1221 


Fig.  19.  Shock  Merging 

These  interesting  properties  of  the  hodographs  reflect  essentially  the  fact  that,  in 
order  to  attain  a  given  particle  velocity,  higher  pressure  has  to  be  reached  by  an  iscntropic 
compression  wave  (equivalent  to  an  infinite  number  of  shocks,  each  propagating  into 
a  regime  produced  by  the  preceding  one)  rather  than  by  a  shock  wave.  In  contrast  to  this, 
the  sound  speed  (or  temperature)  attained  by  the  isentropic  pressure  wave  is  lower  than 
that  produced  by  the  shock  which  produces  the  same  particle  velocity. 

J.4.  Merginc  oi  a  Rarefaclioo  Wave  »ith  a  Shock 

The  various  wave  systems  that  can  be  obtained  by  a  rarefaction  wave  merging  with 
a  shock  have  been  described  by  Glass  and  Hall  (1958).  A  simple  graphic  I  interpretation 
of  these  results  is  given  in  Fig  20  which  illustrates  the  case  of  a  rarefaction  moving  to 
the  right  and  merging  with  a  shock  wave  that  propagates  ahead  of  it.  Referring  to  the 
pressure  velocity  hodograph,  it  appears  that  a  rarefaction  started  from  any  point  of  the 
shock  polar  of  the  same  family  has,  on  one  hand,  a  larger  initial  slope  than  the  polar  but, 
on  the  other,  it  is  straighter  than  the  polar  and,  consequently,  there  must  exist  another 
intersection  such  as  point  2'4'  in  Fig.  20.  Above  this  state,  the  interaction  results  in  a  re- 
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(to  * 

Fig  20.  Merging  of  a  Rarefaction  Wa\c  with  a  Shock 


fleeted  shock,  below  it,  in  a  reflected  rarefaction,  while  state  2'4'  corresponds  to  the  cas^ 
where  there  is  only  a  transmitted  shock  and,  except  for  the  contact  discontinuity,  there 
is  no  reflected  wave.  As  before,  numerical  results  of  the  solutions  arc  listed  in  Table  2. 

4.  Intcractkms  Between,  Shocks  and  Area  Guuiges 

Of  fundamental  importance  to  shock  tube  technology  is  the  interaction  between  a  plane 
shock  front  and  a  section  in  the  tube  where  an  area  change  takes  place.  Basically,  it  is  an 
extension  of  the  classical  problem  concerned  with  the  effect  that  the  change  in  cross-se> 
ction  area  has  on  the  propagation  of  sound  waves  in  tubes.  The  greatest  impetus  to  the 
study  of  this  problem  with  respect  to  finite  amplitude  shock  fronts  was  derived  from 
the  application  of  shock  tubes  as  drivers  for  hypersonic  wind  tunnels  which  involves  the 
interaction  between  the  shock  and  the  divergent  inlet  to  the  test  section.  The  general 
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interest,  however,  into  the  fundamental  features  of  the  interactions  be'  veen  finite  ampli¬ 
tude  waves  and  channels  of  varying  area  led  to  a  number  of  investigations  whirh  included 
other  types  of  area  change,  such  as  convergences  in  the  tube  cross-section,  as  well  as 
obstacles  made  mostly  of  grids  and  wire  gauzes. 

The  main  conclusion  derived  from  these  studies  is  the  realization  that  most  of  the 
experimental  results  check  very  well  with  one-dimensional  analysis  for  a  perfect  gas  with 
constant  specific  heats.  Each  paper  was  concerned  only  with  a  certain  specific  aspect 
of  the  subject  of  shock  interactions  with  area  changes  and,  consequently,  the  coverage 
of  this  field  has  been  so  far  quite  sporadic. 

The  facility  afforded  by  the  vector  polar  method  permits  us  to  attain  a  comprehensive 
view  upon  the  w  hole  subject  of  such  interactions.  The  field  is  covered  by  the  consideration 
of  the  interactions  of  shocks  with,  first,  divergences,  then  convergences  in  the  tube  cross- 
section  area  and,  finally,  with  convergent-divergent  nozzles.  In  the  particular  case  of 
equal  exit  and  inlet,  the  latter  represents  also  an  idealization  of  the  case  of  an  obstacle, 
such  as  a  grid,  affecting  the  whole  cross-section,  that  is  inserted  in  a  constant  area  duct. 

4.1.  Inleraction  witk  Are*  Divergence 

All  the  cases  of  various  plane  wave  systems  that  can  result  from  an  interaction  between 
a  shock  and  an  area  divergence  in  a  tube  are  represented  in  Fig.  21.  Numerical  results 
for  these  problems  are  listed  in  Table  2.  Each  contains  the  traces  of  the  various  wave 
fronts  plotted  to  scale  in  the  non-dimensional  time-space  plane  whose  coordinates  are; 

X  ~  ^  and  T  —  -  where  x  is  the  space  coordinate,  t  is  the  time,  L  is  the  reference 
L  L 

length,  and  Uo  is  the  velocity  of  sound  in  the  undisturbed  medium  which  initially  is  at  rest. 
Shown  underneath  each  wave  diagram  is  a  cut  through  the  wave  system  takon  at  a  suitable 
time  instant.  On  the  right  are  the  graphical  solutions  obtained  by  the  use  of  the  vector 
polar  method  in  the  P  —  L  and  the  A  -  V  planes  for  a  perfect  gas  with  a  specific  heat 
ratio  y  -  1.40. 

Cases  1  through  4  of  Fig.  21  correspond  to  the  same  Mach  number  Mf  =  1.7  of  the 
incident  shock.  The  first  refers  to  a  very  small  divergence  of  an  area  ratio  R  =  1.03.  The 
interaction  under  such  circumstances  produces  eventually  a  rarefaction  wave  traveling 
upstream  and  a  transmitted  shock  which,  as  in  all  the  other  cases  involving  divergences, 
is  weaker  than  the  incident  shock. 

As  the  area  ratio  is  increased,  the  trailing  edge  of  the  rarefaction  fan  is  ’’pushed  back” 
until  it  becomes  stationary  at  the  inlet  to  the  divergence.  Thereupon,  the  interaction  pro¬ 
duces  a  standing  shock  at  some  intermediate  cross-section.  Such  a  situation  is  illustrated 
by  Case  2,  which  corresponds  to  an  area  ratio  R  1.08,  only  a  few  percent  larger  than 
that  of  Case  I. 

When  the  area  ratio  is  increased  still  more,  the  standing  shock  wave  is  ’’pushed  back” 
until  it  finally  becomes  swept  downstream  out  of  the  divergence  producing  the  typical 
configuration  of  hvpersonic  shock  tunnels.  This  is  shown  as  Case  3  corresponding  to 
an  area  ratio  R  1.25. 

Reasoning  ”by  extrapolation”  one  would  expect  at  this  point  that,  as  the  area  ratio 
is  increased  still  more,  the  amount  by  which  the  shock  is  swept  back  would  also  increase 
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indefinitely.  But  this  is  not  the  case  It  turns  out  that  the  amount  of  swecpback  has  an 
upper  bound  so  that,  after  reaching  a  maximum,  it  decreases  until,  with  sufficiently  large 
divergences,  the  snock  Ix'comes  again  stabilized  inside  it  A  u.ise  pattern  obtained  then 
is  demonstrated  by  Case  4  which  corresponds  to  an  area  ratio  R  4.45. 

In  order  to  complete  the  survey  of  the  various  conditions  obtainable  as  a  result  of  the 
interaction  between  a  plane  shock  and  area  divergence,  two  more  solutions  are  shown 
in  Fig.  21  for  the  same  ;  •  'a  ratio  R  -  I  2“)  as  that  of  Case  5,  but  for  different  strengths 
of  the  incident  shock  In  Case  5  the  incident  shock  Mach  number  is  .f/,  1.34,  while 

Case  6  corresponds  to  \f,  }.?.  Conspicuous  in  the  latter  is  the  absence  of  the  rare¬ 

faction  fan  propagating  upstreem  a  charaetcrisiie  feature  of  all  cases  when  the  Mach 
number  of  the  incident  shock  is  high  enough  to  produce  supersonic  flow  in  the  tube.  For 
a  perfect  gas  with  •/  1.4,  such  a  limiting  value  of  the  Mach  number  is  M*  2  069. 

Values  of  vuch  critical  Mach  number  for  other  specific  heat  ratios  are  given  in  Table  L 

4-2.  Interactions  with  Area  (  unvcr^^nces 

In  the  case  of  shock  impingement  on  an  area  convergence  in  a  tube  there  is  an  inte¬ 
resting  peculiarity  in  that  there  exists  a  regime  where  three  different  solutions  satisfy  the 
compatibility  condition  at  the  same  time.  They  correspond  to  the  cases  of  a  reflected 


Fig.  22  ShtKk  Interactions  with  Area  Convergences 
7  \fi  3  20:  F  0.95  8,  A/,  3.20:  F  0.84 

9.  A/j  3  20:  F  0.84  10.  M,  1.34;  F  0.84 
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shock,  a  standing  shock  and  a  transmitted  shock  without  reflection.  This  was  the  subject 
of  one  of  our  earlier  papers  (Oppenheim,  et  al.  (1959))  which  has  been  followed  by  a  number 
of  publications  that  were  concerned  with  a  more  detailed  analysis  of  this  phenomenon, 
namely  those  of  Rudingcr  (1960)  Brower  (1961)  and  Ciolkowski  (1961).  'Moreover  it 
turned  out  that  our  obsersation  about  the  possibilty  of  establishing  more  than  one  wave 
pattern  as  a  result  of  shock  interaction  with  a  converging  duct  received  independently 
an  e.xpcrimental  verification  (Bird  1959). 

for  the  sake  of  completeness,  three  other  representative  cases,  each  for  a  different 
area  ratio  but  all  corresponding  to  the  incident  shock  Mach  number,  A/,  —  3.2.  of  the 
last  example  of  F  ig.  21,  arc  shown  in  Fig  22.  The  first  (Case  7)  for  R  0  95,  represents 
an  interaction  without  reflection,  the  second  (Case  8)  for  R  0  91  illustrates  the  attain¬ 
ment  of  a  standing  shock,  and  the  third  (Case  9)  for  R  —  0.84,  demonstrates  an  interaction 
producing  a  reflected  shock.  With  reference  to  the  latter,  as  the  Mach  number  of  the 
incident  shock  is  decreased,  the  strength  of  the  rarefaction  fan,  3 — 4a,  becomes  weaker 
until  it  disappears  altogether.  This  is  demonstrated  as  Case  10  corresponding  to  the  same 
area  ratio,  R  0.84  as  Case  9  but  to  a  different  incident  shock  Mach  number,  namely 
,U.  -  1.70,  i.e.  the  same  as  that  used  for  Cases  I  to  4.  Again  numerical  results  of  the  so¬ 
lution  are  listed  in  Table  2. 

4.3.  Summary  of  Shock  Interactionn  with  Sincic  Arts  (  hanec* 

As  it  is  apparent  from  the  preceding  discussion,  the  question  as  to  what  kind  of  a  wave 
system  results  from  a  given  interaction  depends  only  on  the  Mach  number  of  the  incident 
shock  and  on  the  area  ratio  that  specifics  the  geometry  of  the  system. 

This  is  summarized  in  Fig.  23  where  all  the  various  regimes  of  solutions  arc  given  in 
terms  of  the  area  ratio  and  the  Mach  number  of  the  incident  shock  for  a  perfect  gas  with 
a  constant  specific  heat  ratio,  y=  1.40.  The  conditions  of  the  various  examples  of  Figs.  21 
and  22  are  indicated  in  Fig.  23  by  points  numbered  correspondingly  to  the  cases  which 
they  represent. 

By  reference  to  Fig.  2.'  one  can  quickly  ascertain  what  kind  of  a  wave  pattern  can  be 
eventually  produced  by  a  given  interaction.  Each  wave  system  occupies  a  regime  in  the 
M,  R  plane.  The  boundaries  have  been  determined  analytically  and  their  properties 
arc  described  in  Appendix  I. 

Fig.  23  can  be  generalized  easily  to  other  values  of  the  specific  heat  ratio,  */  if  the 
effect  of  this  parameter  on  its  characteristic  points  arc  known.  Such  information  is  given 
in  Fig.  24  where  asymptotes  for  lines  III  (Lim  Rm),  IV(Lim  R^y),  l(LimAf,),  and 
V(L  im  \/^)  are  expressed  in  terms  of  •/.  All  these  lines  stem  fre -n  a  point  on  the  ~  1 
axis  representing  the  state  when  the  flow  immediately  behind  the  incident  shock  becomes 
sonic  in  the  tube;  consequently  the  dependence  on  y  of  the  Mach  number  of  such  shocks, 
)/*.  is  also  included  there. 

Since  the  shape  of  the  lines  remains  essentially  the  same,  they  can  be  drawn  for  any 
value  of  '/  on  the  basis  of  the  information  provided  by  Fig.  24.  In  order  to  describe  the 
extent  by  which  '/  can  influence  the  boundaries,  curves  of  Fig.  23  are  redrawn  in  Fig.  25, 
for  two  extreme  values  of  this  parameter,  namely  1.0  and  5/3.  It  is  of  interest  to  note  that 
for  convergences  the  regime  of  ’’peculiarity*’  decreases  quite  rapidly  a< '/  becomes  larger. 
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Fig  23.  Regimes  of  Solutions  for  Sh<Kk  Interactions  with  Single  Area  Changes  in  the  Mach  Niunbcr-Arca 

Ratio  Plane  (7  1.4) 


Fig.  24.  Asymptotes  of  Boundaries  in  the  Area  Ratio-Mach  Number  Plane  and  the  Critical  Shock  Mach 

Number  as  a  Function  of  Specific  Heat  Ratio 


1.-9] 


VKCTOH  POLAh  MKTJIOD 


469 


Fig.  25.  Effect  of  Specific  Heat  Ratio  on  the  Mach  Number-Area  Ratio  Plane,  Fig  2.V 


4.4.  Intcracliont  with  Contt  (anl-DiverKriit  Noiilri 

Several  solutions  of  shock  inte'^actions  with  convergent-divergent  nozzles  are  presented 
in  Fig.  26  and  tabulated  in  Table  2. 

The  first  two  refer  to  the  same  nozzle  whose  geometry  is  fixed  by  area  ratios  /?, 

0.444/1.666,  under  the  influence  of  different  incident  shock  Mach  numbers  1.76 

and  A/,  =  2.581  respectively,  the  former  generating  a  subsonic  and  the  latter  a  supt.  mic 
flow  through  the  duct.  The  resulting  wave  system  of  Case  I  resembles  the  left  side  of 
Case  9  in  Fig.  22  and  the  right  side  of  Case  3  in  Fig.  21.  while  Case  2  in  Fig.  26  is  a  com¬ 
bination  of  the  left  side  of  Case  10  in  Fig.  22  and  of  the  right  side  of  Case  6  in  Fig.  21. 

The  rest  of  solutions  in  Fig.  26  refer  to  experiments  actually  reported  in  the  literature. 
Cases  3  and  4  correspond  to  conditions  produced  by  a  30  mesh  grid  of  38"o  porosity 
taht  has  been  used  by  Franks  (1957)  in  a  shock  tube  w'ith  incident  shock  Mach  numbers 
A/j  1.649  and  2  30  respectively.  Cases  5  and  6  are  those  of  experiments  performed 
by  Dosanjh  (1955)  who  used  a  grid  with  a  solidity  of  25°o  with  incident  shock  Mach  num¬ 
bers  of  A/,  1.58  and  2.30. 

The  results  are  summarized  in  Table  3  which  demonstrates  that,  in  addition  to  its 
simplicity  and  universality,  the  vector  polar  methods  offers  a  sufficiently  good  accuracy 
for  all  practical  purposes. 

Similarly  as  for  single  area  changes,  all  solutions  for  nozzles  can  be  summarized  in 
a  single  A/, — R  plot,  as  shewn  in  Fig.  27  for  y  ^  1.40.  If  the  flow  through  the  convergent 
section  of  the  nozzle  is  wholly  subsonic,  the  resulting  flow  system  can  be  calculated,  di¬ 
sregarding  completely  the  geometry  of  the  nozzle  and  substituting  for  it  only  a  single 
area  change  that  is  specified  by  the  overall  ratio  of  the  exit  to  inlet  cross-section  area. 


Fig.  26.  Shock  Interactions  with  Convergent-Divergent  Nozzles 
I.  A/j  =  1.767:  Rc  =  0.444:  Rj  =  1.666  2.  A/,  =  2.581 :  R^  =  444.0:  R^ 

3.  A/,  =  1.649:  =  0.38:  Rj  =  2.64  4.  A/,  =  2.30:  =  9.38: 

5.  V/j  =  1.58:  Rc  0.75:  =  1.334  6.  Mi  =  2.30:  Rc  =  0.75 
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Such  circumstances  prevail  in  the  regime  to  the  left  of  both  branches  of  Line  I.  To  the  right 
of  this  boundary  the  solution  for  shock  interaction  with  a  no/vlc  can  oc  obtained  by 
combining  solutions  for  consergent  and  divergent  sections.  For  this  purpose  one  has  to 
know  only  the  compatible  slate  at  the  throat  that  represents  the  outflow  conditions 
for  one  and  the  inflow  for  the  other.  It  turns  out  that  such  conditions  are  given  b\  lines  IV 
and  III  in  I  ig.  27.  The  required  matching  of  the  coinergent  part  with  the  discrgence  is 
indicated  there  by  thin  guidelines.  They  connect  corresponding  points  on  the  two  sides 
of  the  R  I  axis  so  that,  having  entered  the  diagram  with  a  given  value  of  .\f,  for  the 
convergence,  the  corresponding  value  of  \f,  for  the  divergence  is  obtained  by  following 
the  appropriate  guideline  up  to  Line  I  or  III  as  illustrated  in  Fig.  27. 


Fig.  27.  Regimes  od  Solutions  for  Shock  Interactions  with  Convergent-Divergent  Nozzles  in  the  Mach 

Number-Area  Ratio  Plane  (y  =  I  40) 

As  previously,  boundaries  of  the  different  regimes  of  solutions  in  the  .\f,  —  R  plane 
have  been  determined  anaivtically  and  arc  described  in  Appendix  2. 

Finally  it  should  be  noted  that  similar  diagrams  can  be  easily  constructed  for  other 
values  of  ;■  by  the  use  of  the  auxiliary  plot  of  Fig.  24  in  the  same  manner  as  it  has  been 
done  for  Fig.  25. 

Incidentally,  the  case  of  a  divergent-convergent  duct  has  been  left  out  of  our  consi¬ 
derations,  since  it  does  not  lead  to  any  interesting  results.  As  contrasted  to  convergence, 
the  divergence  produces  a  change  of  stale  away  from  choking  so  that,  for  most  initial 
conditions,  the  net  results  of  shock  interaction  with  such  a  duct  is  the  same  as  that  pro- 
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duced  b>  a  single  area  change  Any  other  case  can  be  sohed,  of  course,  uiihout  an>  di¬ 
fficulty  by  the  use  of  the  sector  polar  method  follossing  the  same  procedure  as  that  de¬ 
scribed  here  for  the  consergent-divergent  nozzle. 


5.  Performance  of  Shuck  IuIk-s  miIIi  Vrea  CliaiiKC  at  llie  Diaphragm  Section 

One  of  the  well  known  techniques  of  increasing  the  performance  of  shock  tubes  is  the 
use  of  an  area  change  at  the  diaphragm  section  This  has  been  first  proposed  by  1  uka- 
siewicz  (1950).  while  later  Alpher  and  White  (I95!<)  presented  an  analytical  method  for 
the  computation  of  the  performance  of  such  shock  lubes. 

Considered  here  first  arc  cases  insohing  a  single  discrgencc  at  the  diaphragm  section, 
then  a  single  consergence,  and  finally  a  ctrnvcrgent-divcrgenl  nozzle.  As  before,  our  ana¬ 
lysis  leads  to  a  specific  delineation  of  regimes  of  solutions  on  a  chart  with,  this  time,  the 
diaphragm  pressure  ratio  and  the  area  ratio  as  coordinates. 

It  appears  that,  by  increasing  the  diaphragm  pressure  ratio  in  the  case  of  a  gi\cn  di- 
sergence,  one  obtains  first  a  ira  smitted  shock  with  a  reflected  rarefaction,  then  a  standing 
shock  solution  followed  by  a  regime  of  swept-back  shocks,  and  finally,  a  circumstance 
that  somehow  escaped  the  notice  of  others,  of  swept-back  rarefactions.  In  the  case  of 
a  single  convergence  there  arc  only  two  regimes,  one  with  transmitted  shock,  and  the 
other  with  transmitted  shock  and  a  rarefaction  anchored  at  the  inlet  to  the  coinergcnce; 
both  are  associated  with  the  usual  rarefaction,  characteristic  of  the  shock  tube  problem, 
that  propagates  into  the  driver  section.  The  effect  of  the  addition  of  a  convergence  to  pre¬ 
cede  the  divergence,  resulting  in  a  convergent-divergent  nozzle,  is  to  displace  the  boundary 
curves  on  the  diaphragm  pressure-area  ratio  plane  in  the  direction  of  smaller  pressure 
ratios. 


5  I.  l-ffvils  of  Sin|!l('  (  hanevst 

Cases  I  to  4  in  I  ig.  28  are  representative  of  solutions  obtained  with  the  area  divergence. 
The  method  of  present  ition  is  the  same  as  before,  i.e.  each  figure  coiilains  on  the  left 
the  sketch  of  the  physical  A  -T  plane  with  the  wave  world-lines  resulting  in  each  case 
drawn  to  scale,  on  the  right  are  the  hodograph  stale  diagrams  with  the  pressure-velocity 
solution  above  and  the  velocity  of  sound-hodograph  plane  below.  The  substance  in  all 
cases  of  big.  28  is  a  perfect  gas  with  ;■  1.4. 

Case  I  corresponds  to  the  diaphragm  pressure  ratio  /'  2.62  and  area  ratio  (low 

pressure  to  high  pressure  side)  R  I  25.  It  results  in  a  transmitted  shock  fallowed  by 
contact  surface  and  a  rarefaction  propagating  upstream  -  the  typical  wave  pattern  ol 
a  shock  tube  pr.iblem. 

In  Case  2  theie  appears  a  standing  shock  wave  in  the  divergence,  its  e.vistence  is  necessary 
in  order  to  satisfy  the  compatibility  conditions  requiring  the  hodograph  diagram  in  the 
P  -i  plane  to  be  closed.  The  appearance  of  the  standing  shock  is  a'-sociated  with  the 
generation  of  a  new  state  denoted  by  1 .2  so  that  the  flow  in  the  divergence  is  at  first  super¬ 
sonic  from  state  (3)  to  (2)  where  a  change  of  state  takes  place  across  a  shock  from  2  to  1.2 
followed  by  steady  state,  subsonic  flow  from  1.2  to  I. 
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Fig  Is.  VNasc  S> 'lcnl^  Obtained  with  Single  Area  (  lianges  at  the  Diaphragm  Section 
I  Ihl  R  I  2S  2  />,  «  IIS  t.  />,  1 1  |  :  /C  1.25 

4  f\  67  0.  R  I  25  5.  2  62  R  0.444  6  /»,  1 1  I  :  /{  0.444 


■\s  the  pressure  ratio  is  increased  fitrthei.  the  standinu  shock  is  pushed  downstream 
of  tlie  nozzle  so  that  in  Case  3  it  appears  as  a  swept-back  wa\e  2  la. 

I  or  still  h  gher  diaphragm  pressure  ratios  state  3  is  shitted  upwards  in  the  /’  i  plane 
while  its  velocitv  is  fixed  by  the  condition  of  U  1  on  the  —  I  diagram.  Since  the 
change  from  state  3  to  state  2  is  specified  b\  the  area  change,  conditions  are  finalix  achie\ed 
that  cause  state  2  to  be  shifted  to  the  left  side  of  the  incident  shock  polar  on  the  P — C 
plane.  In  order  to  satisfy  the  compatibility  conditions,  the  change  of  state  between  the 
exit  of  the  divergence  and  that  immediately  behind  the  transmitted  shock  involves  then 
the  action  of  a  rarefaction  rather  than  the  swept-back  shock.  An  example  of  such  circum¬ 
stances  is  presented  by  Case  4.  The  boundary  line  between  the  regime  of  solutions  with 
the  swept-back  shock  and  that  with  the  swept-back  rarefaction,  corresponds  to  the  case 
of  only  a  transmitted  shock  and  a  rarefaction  propagating  upstream  with  its  trailing  edge 
anchored  at  the  inlet  of  the  divergence. 
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Case  5  represents  a  typical  result  obtained  with  a  eonsergence,  practically  identical 
with  the  elementary  shock  tube  problem. 

When  the  pressure  ratio  is  increased  while  the  area  ratio  remains  unchanged,  there 
appears  a  rarefaction  propagating  upstream  on  the  low  pressure  side  of  the  convergence 
and  anchored  ,it  its  outlet  cross-section,  as  shown  in  Case  6. 

The  \arious  solutions  described  in  Fig.  28  are  tabulated  in  Table  2  and  arc  represented 
by  the  appropriate  regimes  in  the  diaphragm  pressure-area  ratio  plane.  Fig.  29,  for  the 
case  of  specific  heat  ratio  •/  1.4.  Denoted  there  by  numbers  corresp  inding  to  the  cases 

just  described  are  points  specifying  their  initial  conditions.  The  analytical  determination 


t  t  4  e  I  >e  to  40  to  too 


Fig.  29.  Regimes  of  Solutions  for  '’.Single  .Area  Change  .at  the  Diaphragm  Section  in  the  Pressure  R.atio-Area 

Ratio  Plane  (;'  1,4) 
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Fig.  m.  \symptoie'  of  Doiindanes  m  ihe  Pressure  R;iiio-  \rea  R.itio  Plane  as  a  t  unction  of  Specific 

Meal  R.itio 
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Fig.  31.  Wave  Systems  Obtained  with  a  Convergent-Divergent  Nozzle 
Rc  -  0.444:  =  1.25,  at  the  Diaphragm  Section 

1./’4  =  2.46  2.  P,  =  3.67  3.  P,  H  I  4.  P,  =  36.6 
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of  the  boundary  lines  of  Fig.  29  is  described  in  Appendix  3.  The  chart  gives  a  comprehensive 
survey  of  all  conditions  that  can  be  obtained  with  a  single  area  change  at  the  diaphragm 
section  and  covers  specifically  the  range  of  area  ratios  from  O.I  to  10  and  pressure  ratios 
up  to  150:1.  Of  particular  interest  is  the  regime  characterized  by  point  4  which  somehow 
escaped  the  notice  of  other  investigators. 

Fig.  30  gives  values  of  asymptotes  of  the  various  boundary  lines  that  appear  in  Fig.  29, 
as  a  function  of  the  specific  heats  ratio,  so  that  their  effect  could  be  estimated.  Also  in¬ 
cluded  in  Fig.  30  is  the  effect  of  the  sjsecific  heat  ratio  on  P^,  the  point  of  intersection 
of  all  boundary  lines  in  Fig.  29. 

5.2.  Effect  of  Converitent-Di>erteiit  Noiile 

Fig.  31  represents  the  various  solutions  obtained  for  the  same  convergent-divergent 
nozzle  =  0.444/1.25,  with  different  diaphragm  pressure  ratios.  Numerical  results 

of  these  solutions  are  listed  in  Table  2. 

Case  1  of  /*  =  2.46  gives  the  same  wave  pattern  as  the  elementary  shock  tube  problem. 

Increasing  the  diaphragm  pressure  ratio  to  3.67  leads  to  Case  2  associated  with  the 
appearance  of  a  standing  shock  in  the  divergence.  Similarly  as  in  Case  2  of  Fig.  28, 
it  should  be  noted  that,  as  long  as  choking  is  not  attained,  i.e.  the  line  of  Mach  number 
of  unity  in  the  A  —  V  plane  is  not  reached,  the  solution  is  identical  to  that  described  pre¬ 
viously.  Now  however  the  solution  is  determined  completely  by  the  overall  area  ratio 
(outlet  to  inlet)  of  the  nozzle  rather  than  by  its  shape.  The  convergent  and  divergent  se¬ 
ctions  have  to  be  considered  separately  the  moment  A/  =  1  is  reached  at  the  throat. 


Fig.  32.  Regimes  of  Solutions  for  Convergent-Divergent  Nozzles  at  the  Diaphragm  Section  in  the  Pressure 

Ratio-Area  Ratio  Plane  (y  =  1.4) 
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Case  3  of  I  ig.  31  is  equivalent  to  Case  3  of  Fig.  28  resulting  in  a  sweptback  shock, 
while  Case  4  has  a  similar  wave  pattern  to  Case  4  of  Fig.  28,  namely  the  swept-back  ra¬ 
refaction. 

All  the  results  obtainable  with  a  consergent-clivergent  nozzle  are  summarized  in  Fig.  32 
over  the  same  range  of  the  diaphragm  pressures  and  area  ratios  as  F  ig.  29.  Similar  regimes 
are  obtained  here  as  before,  the  main  effect  of  adding  the  convergence  in  shifting  the 


0  ». 

hig.  .13.  Boundaries  in  the  Pressure  Ratio-Area  Ratio  Plane  for  3  and  !0,  and  the  Critical  Pressure 

Ratio,  P,,  as  a  Function  of  Re- 

boundaries  to  the  left,  that  is  in  the  direction  of  smaller  pressure  ratios.  This  permits 
the  attainment  of  the  same  wave  pattern  with  lower  initial  pressure  ratio  by  the  use  of 
a  nozzle  than  in  the  case  of  a  divergence  alone. 

In  order  to  make  the  coverage  of  the  effect  of  a  nozzle  at  the  diaphragm  more  compre¬ 
hensive,  the  variation  of  the  characteristic  points  of  the  boundary  lines  has  been  computed 
for  the  whole  range  of  possible  convergences  from  0  to  1,  and  the  diaphragm  pressure 
ratios  from  1  to  400  as  shown  on  Fig.  33.  Since  the  shape  of  the  lines  remains  practically 
the  same,  the  det.i  mination  of  the  boundaries  for  any  given  convergence  in  Fig.  32  is, 
with  the  help  of  Fig.  33,  quite  simple. 


5.3.  Case  of  Different  Gaiet  on  the  T»o  Sides  of  the  Diaphrifm 

The  use  of  the  area  change  at  the  diaphragm  section  is  often  associated  with  that  of 
different  gases  on  its  two  sides.  This  has  been  indeed  the  case  treated  by  Alpher  and  White 
(1958)  who  were  concerned  ir  particular  with  helium  and  hydrogen  as  drivers  in  combi¬ 
nation  with  air,  nitrogen  and  argon  in  the  driven  section.  The  problem  then  was  to  eva¬ 
luate,  for  a  given  convergence,  the  pressure  ratio  that  would  drive  a  shock  wave  of  a  given 
strength. 

In  order  to  apply  the  vector  polar  method  to  such  a  problem,  it  is  convenient  to  intro¬ 
duce  an  au.xiliary  diagram.  Since,  as  demonstrated  by  the  graphical  solution  of  Fig.  34 
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>  work  back  from  a  given  slate  1  h>  trial  and  error,  the  effort  is  considerabl\ 
rvdu>.n.l  if  a  single  change  of  state  from  t  to  2  is  known  a  priori.  This  involves  ihe  knowledge 
of  the  change  of  state  produced  by  a  rarefaction  combined  with  a  steady  flow  across  a  duct 
with  a  given  area  change.  Such  a  process  is  lepresented  on  the  auxiliary  graph,  1  ig.  35. 
which  corresponds  to  a  simple  rarefaction  fan  followed  by  steady  flow  through  a  con¬ 
vergence  0.444,  in  the  case  of  the  same  gase.  when  =  I,  and  for  the  helium* 
argon  combination,  corresponding  to  .^4  -  3.16.  The  polars  are  obtained  simply  by 
adding  vecto/ially  the  changes  of  slate  produced  by  a  rarefaction  fan  followed  by  that 
for  a  given  steady  flow  through  the  convergence. 


4 

\ 

\ 


\t  ‘  V  -  W 
0  V  > 


t  - 


« 


Fig.  34.  Wave  System  Obtained  with  a  Convergence  at  the  Diaphragm  Section  LSing  lielium  Driver  and 

Argon  as  Driven  Gas,  P,  70  0;  0.444 


Ik 


\ 

'j-, 

\ 

\ 

f 

t 

xa. 

\ 

\ 

• 

\ 

S  Q4U 

Fig  35.  Auxiliary  Combination  Polars  for  Rarefaction  and  Steady  Flow  throug  Convergence  with 

Rc  0  444  ()•  5  3) 
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A  solution  such  as  that  of  fig,  34  for  I*  70  and  0.444,  is  obtained  with  the  use 
of  the  auxiliary  graph  quite  easily  on  the  /'  U  plane.  In  this  case  the  change  of  state 
from  4  to  2  is  fixed.  For  a  given  position  of  state  I  on  the  initial  shock  polar,  state  2  is 
determined  by  simple  trial  and  error,  so  that  state  4,  at  a  fixed  distance  from  it,  would 
lie  on  the  U  —  0  axis.  Numerical  results  of  the  solutions  are  listed  in  Table  2. 

The  agreement  between  the  solutions  obtained  with  the  use  of  the  vector  polar  method 
.ind  those  obtained  by  Alpher  and  White  is  indeed  quite  satisfactory,  as  demonstrated 
in  Table  3. 


6.  Appendix  I 

Properties  of  Boundaries  in  the  .t/,  R  Plane,  Fig.  23 

6  I  KrKime  of  Kflinied  Rarefaction 


1  inc  I  separates  the  regime  where  the  interaction  produces  a  transmitted  shock  with 
a  rarefaction  fan  prepagating  upstream  and  subsonic  flow  through  the  divergence  from 
that  where  a  standing  sho'  k  appears  inside  the  divergence.  It  is  characterized  by  the  con¬ 
dition  of  sonic  inflow  at  the  throat  with  subsonic  iscntropic  flow  through  the  divergent 
section  The  relationship  between  R  and  3/^  for  this  boundary,  shown  as  Line  1  in  Fig.  23, 
is  determined  from  the  following  equations  describing  the  state  immediately  behind  the 
transmitted  shock  wave  of  Case  I,  Fig.  21: 


(35) 

(36) 


F3 


(>■ 


DC,  2/f, 
r+1 


R 

ty  _  ly 

y-l  y-t  1 


The  velocity  ratio,  L\  —  Ui'ag,  local  velocity  of  sound  .atio,  A^  =  a, -Oo  and  pressure  ratio, 
=  /’,  7’fl.  are,  in  turn,  related  by  the  shock  relations  (12)  and  (13)  while  C,  =  Mg/Oo 
and  /’t  -  P3  arc  also  related  by  equation  (11)  The  Mach  number  of  the  flow  imme¬ 
diately  behind  the  divergence  is  related  to  the  area  ratio  by  the  steady  flow  relation: 


y+i 


which  holds  true  when  the  flow  at  the  inlet  to  the  divergence  is  exactly  sonic  —  the  cha¬ 
racteristic  feature  of  Line  I. 

We  have  thus  seven  equations  for  the  eight  variables  M3  flo-  Ps'Ptx  w,  Mq'  ^1  Pi  Po’ 
\f  ,  and  R,  which  lead  to  a  unique  relation  between  Mi  and  R.  For  /?,  =  x  the 
equations  can  be  reduced  to  a  single  relationship: 

I 

j('/,^  !)(■/  1)  '^  ^  ^  J'(27A//  (r 

(M:  ^\\2yAff 
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which  determines  the  value  of  the  asymptote;  e  g.  for  y  -  1.40  it  yields  Lim  M,  1,15. 

Kl^o. 

Asymptotes  for  other  values  of  •/  are  given  in  Table  I. 

The  incident  shock  Mach  number,  A/^,  at  which  all  boundary  lines  intersect,  is  eva¬ 
luated  from  the  shock  relations  (11),  (12)  and  (13)  for  the  condition  A/,  1  0,  yielding: 


(39) 


4(2—/) 


17)-^ 


1 


which,  for  •/  1.4,  results  in  I/,  2  069.  Critical  Mach  numbers  for  other  values  of  ; 

are  listed  in  Table  I 


6.2.  Regime  of  StAodinc  Shock  Solution 

Line  II  in  F  ig.  23  is  the  upper  bound  for  the  regime  of  the  standing  shock  solution, 
and  represents  the  case  of  a  shock  stabilized  at  the  end  of  the  divergence  F-or  higher 
values  of  \/,  the  flow  velocity  at  the  exit  is  too  high  to  maintain  there  a  stationary  shock, 
and  it  is  therefore  swept  downstream,  as  shown  in  Case  3  of  Fig  21. 

The  boundary  is  determined  by  the  conditions  behind  the  transmitted  shock: 


supplemented  by  the  shock  relations  (I  I),  (12)  and  (13)  giving  the  particle  velocity,  sound 
speed  and  pressure  behind  the  shock  in  terms  of  its  Mach  number.  When  the  Mach  num- 
F>er  of  the  incident  shock  is  lower  than  the  critical  (e  g.  A/,*  -  2.069  for  •/  =  1.4),  the 
flow  into  the  divergence  is  sonic  and  the  relation  between  the  Mach  number  of  the  flow 
immediately  behind  the  divergence  and  the  area  ratio  is  again  given  by  eqn.  (37). 

For  supercritical  Mach  numbers  the  incident  shock  produces  supersonic  flow  with 
respect  to  the  duct  and  the  condition  of  \f  1  at  the  entrance  is  no  longer  valid.  For 
this  reason  equation  (37)  is  not  applicable,  and  the  rarefaction  process  included  in  the 
derivation  of  equations  (40)  and  (41)  is  omitted,  since  there  is  no  rarefaction  wave  pro¬ 
pagating  upstream.  The  new  relations  for  and  arc  then  given  by; 
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Line  II  goes  to  infinit>  in  both  R  and  A/,  so  that  for  an  incident  shock  of  any  strength 
there  is  a  divergence  beyond  \shich  the  interaction  results  in  a  standing  wave. 


6.3.  Regime  of  Peculiarity  of  Shock  Interaction  with  a  Contcrcence 


l  or  small  area  ratios  the  flow  behind  a  sufficiently  strong  incident  shock  remains 
supersonic  across  the  consergence.  Lor  larger  area  ratios  the  same  incident  shock  pro¬ 
duces,  on  the  iMher  hand,  a  reflected  wave  with  subsonic  flow  through  the  convergence 
and  a  choking  state  at  the  throat. 

The  upper  bound  for  supersonic  flow.  Line  III,  in  Lig,  23  corresponds  to  the  case  of 
a  rarefaction  anchored  at  the  exit  of  convergence;  the  lower  bound  for  the  reflected  shock 
solu.ion.  Line  IV,  in  Fig.  23  is  that  of  a  standing  shock  at  the  entrance  to  the  convergence. 
It  turns  out  that  the  two  boundaries  overlap  each  other,  so  that  in  effect  the  R  \f,  plane 
IS  folded  over,  producing  the  regime  of  peculiarity,  referred  to  in  the  text,  where  three 
solutions  are  possible  at  the  same  time. 

The  relationship  for  Line  III  is  simply  obtained  by  combining  equation  (11)  with  the 
equation  for  steady  state  choked  flow  ; 


(44) 

1 

1  2 

f  1 

R  ~ 

uXy  1 

\  2  l\ 

For  Line  IV  the  set  of  equations  is  modified  in  that  the  Mach  number  A/,,  appearing 
in  equation  (44)  instead  of  being  determined  by  the  flow  behind  the  incident  shock,  is 
attained  by  a  normal  shock  propagating  into  this  state,  as  shown  in  Case  9,  Fig.  22.  The 
modification  involves  the  use  of  in  equation  (44)  which  is  related  to  A/,  by  the  normal 
shock  relation; 


(45) 


.\L 


■  «f 

7—1 


The  asymptotic  values  corresponding  to  the  incident  shock  Mach  number  of  infinity 
are  given  by; 


(46) 

(47) 


1 


For  7  —  1.40  this  results  in  Lim  =  0.648  and  Lim  —  0.838  Bo*b  asymptotes 
for  different  values  of  y  are  listed  in  Table  I. 


6.4.  ReKimr  of  Refit'cted  Shock  and  Subsonic  Flow  ihrouch  (onvrrgfiKe 

Line  V  in  Fig.  23  separates  the  regime  where  the  interaction  of  a  weak  shock  with 
an  area  convergence  results  in  a  wave  system  consisting  of  a  transmitted  shock,  a  reflected 
shock  and  fully  subsonic  flow  through  the  convergence  (Case  10,  Fig.  22)  from  that  where 
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chokiny  occurs  at  the  throat  with  an  anchored  rarefaction  there  (Case  9,  f  ip.  22).  It  is 
characterized  by  the  condition  of  sonic  flow  at  the  throat  section  without  the  appearance 
of  the  rarefaction  downstream  of  the  conserpence.  The  boundary  is  obtained  apain  from 
equations  describinp  the  state  behind  the  transmitted  shock  where,  in  this  case,  one  has 
to  take  into  account  the  chanpe  of  stale  broupht  about  by  the  reflected  shock,  so  that  now 


(48) 

(49) 


(/; 


v-i 

(/,  \ 

R  I 


P. 


iy 

V 1 1 


The  \ariables  T,,  /*,,  and  /(,  are  related  to  the  Mach  number  of  the  incident  shock 
wave,  4/,,  b>  equ.itions  (11),  (12),  and  (13)  while  for  sonic  flow  at  the  throat  M,  is  re¬ 
lated  to  R  b\  equation  (43).  Combining  equations  (48)  and  (49)  by  means  of  equation  (II) 
yields  then  one  equation  with  three  unknowns,  M,.  M,,  and  R.  The  additional  relation 
required  to  plot  line  N'  is  supplied  by  writing  i/,  as: 


(50) 


(If  fl, 


which,  by  again  using  the  shock  relations  (II).  (12),  and  (13).  can  be  also  expressed  in 
terms  of  4/,,  and  R.  The  asymptotic  value  of  the  Mach  number  of  the  incident  shock 
wave,  4f,,  which  is  attained  for  the  area  ratio  of  R  ~  0.  is  found  by  imposing  the  addi¬ 
tional  condition  for  the  reflected  shock  from  a  closed  end, 


For  the  case  of  a  perfect  gas  with  y  -  I  40,  the  asymptotic  value  of  4f,  —  1.719.  The 
values  of  Lim  4/,  for  other  values  of  are  shown  in  Table  I 


7.  .Appendix  II 

I'roperlios  of  llotindarirs  in  (he  M,  R  Plane,  hit;.  27 

T  I  Ht'Kime  of  Suhsonu  (  nchoki’d  hinw 

Choking  at  the  throat  of  a  convergent-divergent  nozzle  is  of  particular  importance, 
in  that  it  fixes  the  inlet  condition  for  the  divergent  section  I  or  each  incident  shock  there 
is  a  certain  area  ratio  of  the  convergence  that  produces  choked  flow  without  a  reflected 
shock  wave.  Fine  I  and  \  in  I  ig  27  represent  the  locus  of  such  conditions.  Any  problem 
to  the  left  of  these  lines,  i.e.  when  convergence  does  not  lead  to  choking  at  the  throat, 
may  be  treated  as  a  single  area  disturnance  with  R  R,Rj  representing  the  total  area 
ratio  On  the  other  hand  in  the  regime  to  the  right  of  Fine  V,  there  will  result  choking 
at  the  throat  with  a  reflected  shock,  and  the  corresponding  solution  for  the  divergent 
section  has  to  be  foui.^  by  the  use  of  guide  lines  in  Fig  27. 
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The  equation  for  Line  V  is  deri\ed  b>  combining  the  shock  equations  (II),  and  (12) 
with  the  expression  for  the  area  ratio.  K,  corresponding  to  choked  flow  (44),  yielding 
the  following  closed  form  relation  between  \/,  and  R, . 


(^1) 
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R- 


I) 


/( \f,) 


2 

(I  -  '  n 

r  '  1 
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w  here 


(32) 


/(  '/,) 


(U,-  1)^ 


As  R^  approaches  zero  this  boundary  approaches  its  asymptotic  value  of  unity,  ie  : 


Lim  t/,  1 

-H 

The  equation  for  Line  I  is  obtained  by  noting  that  the  flow  through  the  no/zle  is  isen- 
tropic.  Therefore  the  boundary  l.inc  I  corresponds  to  the  case  of  nozzles  l>ctween  tubes 
of  equal  cross-sectional  area,  i.e  R  R,Rj  I,  and  Line  1  is  simply  a  reciprocal  of 
Line  V.  Hence,  from  inspection  of  equation  (51)  its  equation  is: 


(53)  R-j 

and  the  asymptote  is: 


Lim  M 


7.2.  Recimv  of  Standing  Shoik  Solution  is  the  Onergence 


Beyond  Boundary  Line  I  a  standing  shock  is  generated  somewhere  inside  the  diver¬ 
gence  section  of  the  nozzle,  l  ine  II  serves  the  same  purpose  as  Line  II  in  Lig.  23,  namely 
It  represents  the  case  when  the  standing  shock  created  in  the  divergence  is  stabilized  at 
its  end.  One  can  deduce  the  relation  between  R  and  .V/j  for  Line  II  by  noting  that  once 
the  flow  becomes  choked  at  the  entrance  to  the  divergence,  the  flow  through  it  is  inde¬ 
pendent  of  the  processes  producing  the  choked  condition  To  each  state  abng  the  M  1 
locus  there  corresponds  a  value  of  area  ratio,  Rj.  that  will  produce  a  standing  shock  at 
the  exit.  In  the  case  of  single  area  disturbance  described  in  I  ig.  23  such  state  was  deter¬ 
mined  by  the  change  of  state  produced  by  the  incident  shock,  \/,,  followed  by  a  reflected 
rarefaction  wave.  In  the  nozzle  case  it  corresponds  to  a  different  incident  shock  \f,',  fol¬ 
lowed  by  a  steady  flow  process  through  the  convergence,  lor  the  same  flow  velocity  at 
the  throat,  i  e.  i ~  I  n  the  two  incident  shock  Mach  numbers  are  related  as  follows: 


r-i 


where,  by  substitution  of  equations  (II),  (12),  and  (53),  the  left  hand  side  of  the  equation 
may  be  then  expressed  in  terms  of  Atl  and  the  right  hand  side  in  terms  of  \//.  I  sing  then 
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the  information  of  Fig.  23  on  the  dependence  of  the  area  ratio  for  a  single  divergence 
on  the  incident  shock  Mach  number,  together  with  equation  (54),  one  can  plot  Line  II 
in  Fig.  27. 

For  stronger  incident  shocks,  i.e.  A/,  >  A/^,  the  flow  through  the  convergence 
becomes  supersonic  and  the  compatibility  condition  for  the  divergent  section  dominates 
the  establishment  of  the  entire  flow  process.  The  expression  for  Line  II  can  be  then  deduced 
only  from  the  shock  relations  (11)  and  (12),  similarly  as  it  was  done  for  Line  II  in  Fig.  23. 
The  flow  velocity  and  pressure  ratio  immediately  behind  the  transmitted  shock: 


y-i  )— 1 


are  related  by  equation  (11)  while  /<,,  R^,  and  Afi  can  be  expressed  in  term  of  A/, 

by  means  of  equations  (11),  (12),  (13)  and  (44)  and  Ms  in  terms  of  Rj  by  means  of  equ¬ 
ation  (37).  Consequently,  the  above  set  of  equations  fixes  the  relationship  between 
and  A/j  for  Line  II. 

In  the  limit  when  A/,  -»  x,  the  results  of  Fig.  23  still  applies,  i.e.  for  Line  11. 

Lim  R.^  =  CO 

7.3.  Resime  of  Transmitted  Rarefaction  Propagation  t'pstream 

In  cases  of  small  area  divergences,  stronger  incident  shocks  produce  weaker  auxiliary 
shocks  until  finally  the  latter  are  annihilated  forming  a  rarefaction  fan  ins'tead.  This  tran¬ 
sition  of  shock  wave  into  rarefaction  fan  occurs  along  Line  III.  It  is  determined  from  the 
equation  of  supersonic  choked  flow  through  the  convergence  without  reflection,  followed 
by  supersonic  expansion  through  the  divergence  to  the  state  behind  the  transmitted  wave. 
Line  Ml  is  then  simply  the  isentropic  supersonic  extension  of  Line  1  and  the  relation  bet¬ 
ween  and  Ml  is  again  given  by  equation  (53),  where  A/j  is  now  restricted  to  values 
greater  than  Af^.  As  the  incident  shock  Mach  number  approaches  infinity,  Rj  approaches 
its  asymptotic  value  of 

_ 1  _ 

(57)  Lim  /?,,,.  =  («/-!) 

Mi -00  \2l 

which  for  y  =  1.40  results  in  Lim  Rj  =  1 .543.  The  result  for  other  values  of  y  are  given 

Mj-oo  *" 

in  Table  I. 


7.4.  Regioief  of  Solutions  for  Cobrergent  Soctioo 

Strong  incident  shock  waves  impinging  on  nozzles  with  small  convergences  may  pro¬ 
duce  standing  shocks  inside  the  convergent  section  in  the  same  way  as  for  a  single  con¬ 
vergence.  The  regime  of  such  cases  is  delineated  by  Lines  III  and  IV  identical  to  those  of 
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I  ig.  23.  The  rest  oi'  the  wave  svstem  depends  on  the  divergent  nart  ol  the  1107/le  The 
regime  between  l  ines  l\  and  III  represents  the  case  of  siipersonic-unehoked  flow  at  the 
throat  Similarl)  as  the  regime  between  Lines  V  and  I.  the  resulting  wave  svstem  is  com¬ 
puted  here  again  as  if  there  was  onl>  a  ^ingic  area  disturbance  with  an  area  ratio  equal 
to  the  overall  area  ratio  of  the  nozzle,  1  c  R  R^Rj- 

N.  Appi-ndix  III 

I’ropiTties  <if  liouixLtrii's  in  Ihi'  li  /’  Plane,  I'jg.  29 

N  I  Ki'Kimp  ol  thi'  RolKtU-d  Uurvf n* tion  Uatp  «ith  Subsonu  Hom  throu(;h  ihe  Ditprecnt'p 

Line  I  separates  the  regime  representing  solutions  with  a  tran  mitted  shock,  a  rare- 
factiem  fan  propagating  upstream,  and  subsonic  flow  through  the  divergence,  from  that 
where  a  standing  shivck  wave  appears  inside  the  divergence  It  is  characterized  bv  the 
condition  of  sonic  inflow  at  the  throat  with  iscntropii.  flow  through  the  divergent  section. 
The  expression  between  R  and  /’,  for  this  boundarv,  shown  as  Line  I  in  I  ig.  29,  is  deter¬ 
mined  from  the  following  equations  describing  the  state  immediately  behind  the  trans¬ 
mitted  shock  wave  of  Case  I,  Fig.  28; 

(58)  r,  “  '  ' 


The  Mach  number  of  the  flow  immediately  downstream  of  the  divergence  is  related  to  the 
area  ratio  by  the  steady  flow  expression  (37)  while  equations  (58)  and  (59)  can  be  com¬ 
bined  by  means  of  the  shock  equation  (II),  leading  thus  to  a  unique  relation  between 
R  and  f\. 

The  asymptotic  value  c'f  /’,  corresponding  to  the  area  ratio  of  infinity  is  given  by: 


l.imt/’,)  I  ‘  I 


which,  for  ;■  1.4,  results  in  l.im  (Z’,),  1.89.  (Lor  other  values  of-/  see  Table  I). 

The  diaphragm  pressure  ratio,  /’*,  at  which  all  boundary  lines  intersect,  is  determined 
by  setting  the  are;i  ratio  and  the  Mach  number  of  the  flow  behind  the  incident  shock 
wave  equal  to  1,0,  Lquations  (58)  and  (59)  evaluated  with  the  use  of  equation  (II),  lead  to: 


which,  for  /  1.4,  yields  10.35.  (For  other  values  of  /  sec  Table  I). 


1.2.  Rrcimf  of  Standing  Shoik  Solution 


Line  II  in  Fig,  29  is  the  uppter  bound  for  the  regime  of  the  standing  shock  solution 
and  represents  the  case  of  a  shock  stabilized  at  the  end  of  the  divergence. 


4S6 


K  N  MFIM.  P,  a  I  HTIKW  AND  A.  .1  I. ADI  HMAN 


1461 


The  boiiiul.n)  is  (.lelcrmincd  from  the  conditions  behind  the  transmitted  shock  \\a\e 
of  (  .l^c  V  1  It?  2s : 


uhich,  loiicther  uith  the  shock  equation  (II)  and  the  steady  flow  relation  (37)  between 
the  are.i  ratio  and  the  Mach  number  of  the  flow  inside  the  exit  plane  of  the  divergence, 
vields  the  desired  relation  between  R  and  /\ 

No  asymptotic  value  exists  for  Line  II  since  both  R  and  /’i  go  to  infinity. 

1.3.  Reeimf  of  Tran^mittod  Rarefaction  Propaxating  I  p%tream 

\s  the  standing  shock  wave  at  the  exit  of  the  divergence  is  swept  downstream,  it  de¬ 
creases  in  strength  and  eventually  becomes  a  downstream  swept  rarclaction  wave.  Line  III 
is  the  locus  of  states  representing  the  transition  from  compression  to  rarefaction.  This 
case  is  identical  to  the  wave  system  corresponding  to  Line  I,  with  the  exception  that  the 
flow  behind  the  transmitted  shock  is  now  .supersonic,  so  that  equations  (5S)  and  (59) 
can  be  also  used  to  determine  Line  III.  Similar  to  Line  II.  no  asymptotic  value  exists 
since  again  both  R  and  go  to  infinity. 

8.4.  Regime  of  \nchored  Rarefaction  Downstream  if  I'oniergence 

Line  IV,  Fig.  29,  separates  the  regime  representing  solutions  with  a  transmitted  shock, 
a  rarefaction  fan  propagating  upstream,  and  subsonic  flow  through  the  convergence, 
from  that  where  a  rarefaction  wave  is  anchored  at  the  throat  of  the  convergent  section. 
It  is  characterized  then  by  the  condition  of  choked  flow  at  the  throat  section. 

The  boundary  is  determined  by  the  conditions  behind  the  transmitted  shock  wave 
of  Case  5,  f  ii!.  28: 


y — 1 


which,  together  with  the  Nhock  wave  equation  (II)  and  the  steady  flow  relation  (44)  be¬ 
tween  the  area  ratio  and  the  Mach  number  of  the  flow  into  the  convergence,  arc  sufficient 
to  relate  uniquely  R  to 


VF.CTOH  POhAR  MKTHOD 


487 


H7) 


The  asymptotic  \alue  of  l\  corresponding  to  /?  0  is  given  by: 


(65)  Lim  (Ti),,  - 

R-O 

uhlch,  for  y  1.4,  yields  Li  n  {P^)iy. 

R-O 

of  •/  are  listed  in  Table  1 


5.92.  Asymptotic  values  of  for  other  values 


9.  Appendix  IV 

I’ropertics  of  Boundar)  Lines  in  the  R-P  Plane,  Fig.  32 

The  procedure  for  the  derivation  of  equations  for  the  boundary  lines  of  Fig.  32  is 
essentially  the  same  as  that  used  for  the  case  of  a  single  aica  change  at  the  diaphragm 
section.  Fig.  29  However,  for  the  present  ease  the  diaphragm  pressure  ratio.  P^.  is  a  fun- 
tion  of  the  two  independent  variables  R  and  Rj  which  specify  the  geometry  of  the  noz/lc. 
It  is  necessary  then  to  choose  one  of  the  area  ratios  as  the  independent  variable  and  the 
other  as  a  par.imeter  for  which  a  family  of  boundary  lines  could  be  determined.  As  ap¬ 
parent  from  Fig.  29,  it  is  of  greater  interest  to  consider  R^  as  the  variable  ;ind  R^  as  the 
para  meter. 

Uoundary  lines  for  a  convergent-divergent  noz/le  have  been  computed  for  three  values 
of  R^:  0,  0.444,  and  1.0.  Since,  for  each  boundary  line,  sonic  flow  exists  at  the  throat,  the 
effect  of  the  convergence  is  taken  into  account  by  the  following  relations  describing  the 
flow  velocity  and  local  pressure  at  the  throat  section: 


1 


and  the  derivation  of  the  boundary  lines  is  then  similar  to  that  used  for  the  single  diver¬ 
gence.  For  sonic  flow  at  the  throat,  the  Mach  number,  t/j.  of  the  flow  into  the  convergence 
IS  related  to  the  area  ratio,  R^,  by  equation  (44). 

9.1.  Recime  of  Subtionic  Flow  throuKh  the  *)i«ergeiire 

Line  1  of  Fig.  32  serves  the  same  purpose  as  Line  I  in  Fig.  29,  in  that  it  represents  the 
case  of  sonic  inflow  at  the  throat  with  subsonic  isentropic  flow  through  the  divergence, 
separating  regimes  of  Case  1  from  that  of  Case  2  of  Fig.  31.  In  the  present  instance,  ho¬ 
wever,  the  conditions  behind  the  transmitted  shock  wave,  related  by  the  shock  equation  (1 1), 
are  given  by: 

(68)  L\  -  ’ 

K 

(69)  P,  P, 

Pi 
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where  A/,„  the  Mach  number  of  the  flow  at  the  outlet  of  the  divergence,  is  related  to  Rj 
by  steady  flow  expression  (37). 

I  2  \  ‘ 

Since  -  j  —  -)*<<•-•)  and  =  I  as  Rj  goes  to  infinity,  the  expression  for 

\y  +  1 ' 

the  asymptotic  value  of  can  be  obtained  directly  from  equation  (69) 


Lim  l\ 

R<(-«  pr 


1 


y 

7-1 


where  />4  Pr,  determined  b>  the  area  ratio  of  the  convergence,  is  given  by  equations  (66) 
and  (67).  For  R^  =  1.0  the  asymptotic  expression  reduces  to  that  obtained  for  Line  1, 
Fig.  29,  namely 

Lim  if\)  =  I 

\y  ,  \J 

which,  for  y  ^  1-40,  yields  Lim  /’,  1.89,  while  for  R^  —  0, 

Fj  —ru 

Lim  ^4  =  1  0 

-•'X 


9.2.  Rfcjm«  of  Standing  Shock  .Solution 

Similarly  as  Line  II  in  Fig.  29,  Line  11  in  Fig.  32  is  the  upper  bound  for  the  regime  of 
the  standing  shock  solution  and  represents  the  case  of  a  shock  stabilized  at  the  outlet  of 
the  divergence  The  boundary  is  determine  J  from  the  shock  equation  (II)  and  the  follow  ing 
expressions  for  conditions  behind  the  transmitted  shock  wave  of  Case  3.  Fig.  31: 


where  A/.,  the  supersonic  critical  Mach  number  at  the  exit  of  the  divergence,  is  related 
to  Rj  by  equation  (37). 

No  asymptotic  value  exists  for  Line  II  since,  as  in  the  case  of  the  divergence  at  the 
diaphragm  section,  both  Rj  and  go  to  infinity, 

9.3.  R  'Kimr  of  TraotiniKrd  Rarefaction  PropaKaiinit  I  pt'  ram 

As  the  pressure  ratio  P^  increases,  the  strength  of  the  downstream  swept  shock  wave 
of  Case  3,  Fig,  32,  decreases  and  eventually  becomes  a  downstream  swept  rarefaction 
wave.  Line  III,  similar  to  Line  III  of  Fig.  29,  is  then  the  locus  of  states  representing  the 
transition  from  compression  to  rarefaction.  The  wave  system  for  Line  111  is  therefore 
identical  to  that  for  Line  I  with  the  difference  that  now  supersonic  flow  exists  at  the  exit 
of  the  divergence.  Equations  (68)  and  (69)  can  be  also  used  to  determine  Line  111,  except 
that  the  critical  Mach  number  at  the  end  of  the  divergent  section  is  evaluated  for  super¬ 
sonic  flow. 

Similarly  as  for  Line  II,  no  asymptotic  value  exists  for  Line  III,  since  again  both  Rj 
and  P^  go  to  infinity. 
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Table  2 

Solutions  of  Wave  Interactin  Problems 


1 

State 

M 

P, 

An 

t/. 

y 

H' 

-1 

Fig.  17 

4 

— 

2.62 

1 

_ 

* — 

1 

la 

-1 

1.584 

0.933 

0.33 

-1 

-1 

-1 

-  0.933 

-0.60 

1 

1.25 

1.584 

1.067 

0.33 

1.25 

1.25 

2 

-1.2 

2.41 

1.133 

0 

-1.28 

-  0.95 

3 

•  1 

0.887 

0  86 

0 

0.933 

1.263 

+  1 

0.86 

0.86 

Fig.  18 

1 

1.6 

2.75 

1.173 

0.8 

1.6 

1.6 

10 

-0.14 

0.61 

2.47 

1.8 

-0.14 

-  1.4 

2 

-0.08 

1.336 

2.44 

3.02 

-0.21 

0.59 

20 

-1 

1.15 

2.41 

3.28 

-0.244 

0  58 

-1 

-0.241 

0.87 

20a 

1.37 

1.15 

2.62 

3.28 

3.38 

5.18 

Fig.  19 

1 

1.6 

2.77 

1.176 

0.80 

1.6 

1.6 

2 

1.6 

7.7 

1.39 

1.77 

1.88 

2.68 

3a 

-1 

7.24 

1.377 

1.83 

-1.39 

0.38 

1 

( 

-1 

-1.377 

0.453 

! 

3 

2.54 

7.24 

1.49 

1.83 

2.54 

2.54 

Fig.  20 

1 

3.92 

18 

2 

3.1 

3.92 

3.92 

2 

-1 

8.6 

1.8 

2.1 

2.0 

5.1 

-  1 

1.8 

3.9 

3 

-1.02 

.05 

1.81 

2.04 

-1.834 

.27 

4 

2.77 

8.9 

1.57 

2.04 

2.77 

2.77 

2' 

1 

3.56 

1.585 

1.04 

2.0 

5-'  1 

-1 

1.585 

2.625! 

r 

— 

— 

_ 

— 

_ 

4' 

1.78 

3.56 

1.246 

1.04 

1.78 

1.78 

2" 

1 

1.65 

1  425 

0.2 

2.0 

5.1 

•  1 

1.425 

1.625 

. 

1 

3" 

-1 

1.506 

1.405 

0.32 

- 1 .425 

-  1.225 

-1 

-1.405 

-1.085 

4" 

1.2 

1.506 

1.067 

0.32 

1.2 

1.2 

Fig.  21—1 

1 

1.7 

3.21 

1.207 

0.927 

1.7 

1.7 

R  =  1.03 

2 

-1 

2.72 

1.173 

1.08 

-1.207 

-0.280 

I  * 

-1 

-1.173 

-0.093 

il 

3a 

— 

3.06 

1.201 

0.89 

— 

— 

1.668 

3.06 

1.201 

0.89 

1.668 

1.668 

1 

1.7 

3.21 

1.207 

0.927 

1.7 

1.7 

-1 

-1.207 

-0.280' 

Fig.  21—2 

2 

-1 

2.45 

1.161 

1.161 

-1  161 

®  1 

R  =  1.08 

3 

— 

1.792 

1.11 

1.393 

— 

:  A/(  =  1.7 

3.4 

-1.25 

2.97 

1.195 

0.977 

-1.388 

0.005 

4a 

— 

3.06 

1.203 

0.907 

— 

— 

1 

4 

1.683 

3.06 

l.a)3 

0  907 

1.683 

1.683 

4* 
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Table  2  (continued) 


State 

A/ 

Po 

A, 

1/. 

y 

tv  ! 

1 

1.7 

3.21 

1.207 

0.927 

1.7 

1.7  ' 

-1 

-1.207 

-0.280! 

Fig.  21—3 

2 

-1 

2.45 

1.169 

1.161 

-1.161 

0 

R  =  1.25 

3 

— 

1.093 

1.033 

1.671 

— 

— 

Af,  =  1.7 

4a 

-1.578 

2.99 

1.214 

0.87 

-1.63 

0.0411 

4 

1.65 

2.99 

1.197 

0.87 

1.65 

1.65 

1 

1.7 

3.21 

1.207 

0,927 

1.7 

1.7 

-1 

-1.207 

-0.280 

Fig.  21—4 

2 

-1 

2.45 

1.169 

I.I6I 

-1.161 

0  ! 

R  =  4.45 

3 

— 

3.21 

0.765 

2.20 

— 

— 

.V/i  =  1.7 

3.4 

-2.876 

1.786 

1.285 

0.60 

-2.20 

0 

4a 

— 

1.878 

1.244 

0.47 

— 

— 

4 

1.320 

1.878 

1.0975 

0.47 

1.320 

1.320 

1 

1..U 

1.955 

1.105 

0.5 

1..34 

1..14 

-1 

-1.105 

-0.605 

Fig.  21—5 

2 

-1 

1.70 

1.079 

0.64 

- 1 .079 

-0.439 

R  =  1.25 

3a 

— 

1.896 

1.096 

0.46 

— 

— 

Ml  =  I..34 

1.33 

1.896 

1.096 

0.46 

1.33 

1.33 

1 

3.2 

11.7 

1.71 

2.40 

3.2 

1 

3.2  ! 

Fig.  21 — 6 

3 

6.83 

1.583 

2.80 

— 

— 

R  =  1.25 

4a 

-1.219 

10.7 

1.70 

2.26 

-1.93 

0.87 

Mt  =  3.2 

4 

2.%2 

10.7 

1.659 

2.26 

2.962 

2.962 

Fig.  22-7 

1 

3.2 

11.7 

1.71 

2.40 

3.2 

3.2 

R  =  0.95 

3 

14.1 

1.745 

2.23 

— 

— 

-1 

-1.745 

0.485 

Mi  =  3.2 

4a 

-1 

12.2 

1.714 

2.44 

-1.714 

0.726 

4 

3.26 

12.2 

1.726 

2.44 

3.26 

3.26 

1 

3.2 

11.7 

1.71 

2.4 

3.2 

3.2 

Fig.  22—8 

1.2 

— 

14.1 

1.745 

2.23 

— 

— 

R  =  0.91 

2 

-1.28 

24.3 

1.90 

1.52 

-2.23 

0 

•II 

3 

— 

19.3 

1.845 

1.845 

— 

— 

-1 

-1.845 

0 

4a 

-1 

12.4 

1.722 

2.46 

-1.722 

0.738 

4 

3.279.12.4 

i 

1.738 

2.46 

3.279 

3.279 

1 

11.7 

1.71 

2.40 

3.20 

3.20 

Fig.  22  9 

2 

-1.48 

28.5 

1.962 

1.20 

-2.53 

-0.13  , 

R  =  0.84 

3 

— 

19.3 

1.87 

1.87 

_ 

_  ! 

-1 

-1.87 

0  i 

Mi  =  3.2 

4a 

-• 

12.55 

1.74 

2.48 

-1.74 

-0.74 

4 

3.30 

12.55 

1.746 

2.48 

3.30 

3.30 

Fig.  22—10 

1 

1.70 

3.21 

1.207 

0.927 

1.70 

1.70 

/I  =  0.84 

2 

-1.121 

4.18 

1.249 

0.70 

-1.354 

-0.427 

Mi  =  1.7 

3a 

— 

3.35 

1  208 

0.97 

— 

3 

1.737 

3.35 

1.221 

0.97 

1.737 

1.737 
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Table  2  (continued) 


State 

M 

Aq 

f/o 

V 

1 

1.767 

3.44 

1.228 

1.  3 

1.767 

1.767i 

Fig.  26—1 

2 

-1  373 

7.1 

1 .356 

0.36 

-1.688 

-0.688 

r 

— 

3.94 

1.24 

1.24 

— 

—  1 

Ri  =  1.66 

3 

— 

0.97 

1.017 

1.96 

— 

— 

Mi  =  1.767 

4a 

-  1.754 

3.32 

1.258 

0.96 

-1.785 

0.175 

1 

4 

1.729 

3.32 

1.215 

0.96 

1.729 

1.729 

: 

' 

1 

2.581 

7,60 

1.489 

1.83 

2.581 

1 

2.581 

'  Fig.  26.2 

2 

-1.681  23.9 

1.79 

0.50 

-2.50 

-0.67 

Rc  =  OM 

T 

— 

13.3 

1.635 

1.63 

— 

1 

Ri  =  1.66 

3 

_ 

3.285 

1.339 

2.63 

— 

1 

Mi  =  2.58 

4a 

-1.431 

7.3 

1.53 

1.78 

-1.915 

0.715, 

i 

4 

2.530 

7.3 

1.47 

1.78 

2.53 

2.53  ! 

1 

1 

1.649 

3.00 

1.192 

0.868 

1.649 

1.649 

Fig.  26—3 

2 

-1.304 

5.48 

1.304 

0.30 

-1.56 

-0.692 

!  Rc  =  0M 

T 

_ 

2.98 

1.20 

1.20 

— 

Ri  =  2.64 

3 

— 

0.331 

0.87 

2.07 

— 

_  t 

1  Mi  =  1.649 

4a 

-2.419 

2.18 

1.231 

0.60 

-2.1 

-0.03  ‘ 

1 

4 

1.418 

2.18 

1.125 

0.60 

1.418 

1.418 

1 

1 

2  30 

6.00 

1.395 

1.554 

2.30 

! 

2.30  1 

Fig.  26 — 4 

2 

-1.61017.1 

1.642 

0.37 

-2.241 

-0.6871 

!  J?c  =  0.38 

r 

9.3 

1.506 

1.506 

_ 

1 

;  Ri  =  2.64 

3 

1  02 

1.094 

2.78 

!  Mi  =  2.30 

4a 

-2.063 

4.9 

1.445 

1.34 

-2.26 

0.52  ! 

i 

1 

4 

2.084 

4.9 

1.325 

1.34 

2.084 

2.084 

1 

■ 

1 

1 

1.58 

1 

2.745  !  1.172 

0.789 

1.S8 

1.58  ! 

Fig.  26—5 

2 

-1.105 

3.46 

1  209 

0.59 

-1.295 

-0.50 

1  Rc  =  0.75 

r 

— 

2.18 

1.131 

1.12 

— 

—  I 

1  Ri  =  1.334 

3 

— 

0.827 

0.986 

1.70 

— 

—  j 

'  Mi  =  1.58 

4a 

-1.716 

2.66 

1.20 

0.77 

-1.695 

0.005 

1 

4 

1.557 

2.66 

1.166 

0.77 

1.557 

1.557; 

1 

2.  .30 

6.00 

1.395 

1.554 

2.3 

2.3 

Fig.  26-  6 

2 

-1.393  12.6 

1.555 

0.78 

-1.945 

-0.391 

i?c  =  0.75 

T 

— 

791 

1.454 

1.44 

-  1 

Ri=  1.334 

3 

— 

3.00 

1.267 

2.14 

— 

Mi  =  2.3 

4a 

-1.348 

5.86 

1.411 

1.52 

-1.708 

0.432' 

1 

4 

2.273 

5.86 

1.375 

1.52 

2.273 

2.273 

Fig.  28  -1 

4 

- 

2.62 

_ 

i 

—  1 

-1 

-1.0 

-1.0 

R  =  1.25 

3 

-1 

1.448 

0.917 

0.41 

-0.917 

-0.507: 

P  =  2.62 

la 

— 

1.537 

0.926 

0.30 

— 

1 

1 

1.208 

1.537 

1.06 

0-0 

1.208 

1.208, 
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Table  2  (continued) 


State 

A/ 

_ 1 

P. 

.4o 

V 

w 

4 

_  1 

546 

_ 

— 

-1.0 

-I.O 

Fig.  28-2 

3 

— 

1.537 

0.8 

0.83 

-0.83 

0 

R  1  25 

2 

— 

0.772 

075 

1  14 

- 

— 

P  ---  5  47 

1.2 

-1.52 

1  952 

0  87 

0  603 

-1.14 

0 

la 

2.042 

0  88 

0.539 

— 

1 

1.373 

2042 

1.113 

0  539 

1.373 

1.373 

4 

II. 1 

— 

— 

-10 

-1.0 

Fig.  28-3 

3 

—  1 

3  125 

0.83 

0.83 

-0.83 

0 

R  =  \. 25 

2 

— 

1  377 

0  741 

1.196 

— 

— 

P  11.1 

la 

-1.353 

2.72 

0823 

0  78 

-1  00 

0.196 

1 

-1.576 

2  72 

1.173 

0.78 

1.576 

1.576 

4 

_ 

67  0 

— 

— 

— 

— . 

-  1 

-1.0 

-1.0 

Fig  28  -4 

3 

—  1 

18.7 

0  83 

0.83 

083 

0 

«  -  1.25 

2 

— 

8.25 

0.741 

1.196 

— 

— 

—  1 

-0.741 

0.455 

P  -  61 

1  la 

—  1 

5.31 

0.693 

1.42 

-0  693 

0.727 

1 

2  166 

5.31 

1.352 

1.42 

2.166 

2.166 

4 

_ 

2.62 

— 

_ 

—  1 

-1.0 

-1.0 

Fig.  28-5 

3 

—  1 

2.115 

096 

0  20 

-0  96 

-0  76 

R  0.444 

la 

— 

1  822 

0  943 

0  46 

P  2.62 

1 

1  312 

1,822 

1  096 

046 

1.312 

1.312 

4 

III 

- 

- 

— 

- 1 

-1.0 

10 

Fig.  28  6 

3 

—  1 

8  4 

0951 

0.25 

-095 

0  70 

R  0444 

2 

— 

4  71 

0  883 

0  883 

—  1 

-0.883 

0 

P  11  1 

la 

1 

3  74 

0  848 

1.08 

-0848 

0.232 

1 

1  8.3 

.3.74 

1.252 

1  08 

1  83 

1.83 

Fig.  }\-  1 

4 

2.46 

— 

— 

R^  0.444 

la 

-  1 

1.70 

0  944 

040 

1 

Rj  1.25 

1 

1.265 

1.70 

0925 

0  40 

1  265 

1.26'i 

P  2.46 

4 

1 

1 

3.67 

_ 

1 

Fig.  31  2 

3 

1 

1 

2  51 

0  947 

0,24 

1 

-0  947 

—  1 

-0.707 

R,  0.444 

- 

1.405 

0.87 

0.807 

— 

— 

Ri  1  25 

2 

— 

0.764 

0.795 

1  16 

— 

— 

P  3  67 

1  2 

1.46 

l.8(M 

0.914 

0.63 

-1.16 

0 

la 

2  075 

0  932 

055 

- 

— 

1 

-1.38 

2075 

1  116 

0  55 

1.38 

1  38 
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Table  2  (continued) 
State  M  '  t/,  y  If'  ' 


4 

II 

1 

- 

1 

-- 

-1 

-0.707 

Fig.  .^1-3 

3 

—  1 

i  7.7 

0.947 

0  24 

1 

0  947 

=  -0.444 

—  4.22 

0.87 

0  807 

—  1 

Ri  1 .25 

2 

—  1.86 

0.774 

1  24 

- 

R  II. 1 

la 

-1289  3.22 

0.838 

0.93 

0.995 

0.245 

1 

1 

1.704  3  22 

0.83 

0.93 

1.704 

1.704 

4 

36.7 

_ 

_ 

— 

1 

-10 

-10 

Fig.  31  4 

3 

-1  25.00 

0  947 

0  24 

0947 

-0.707 

Rc  0.444 

1400 

0.87 

0.807 

— 

Ri  125 

2 

—  6.18 

0.774 

1  24 

- 

1 

-0.774 

0466 

1  P  36.6 

la 

-1  5  10 

0  748 

I..37 

-0.748 

0.622 

1 

2.124  5  10 

075 

1  37 

2  124 

2.124 

4 

70.0 

3.16 

— 

— 

— 

—  1 

1 

1 

,'.I6 

-3.16 

Fig.  34 

3 

- 1  47.00 

'  2  90 

0.72 

-2.90 

-2  18 

R  =  0.444 

2 

—  24  50 

2.55 

2  55 

— 

—  1 

2.55 

0 

■X! 

11 

o 

la 

-1  20.5 

2.44 

2.88 

-2.44 

0.44 

1 

1 

4.075  20.5 

2.47 

2.88 

4.075 

4.075 

Tab  1 

e  3 

Solutions  Obtained  by  Means  of  the  Vector  Polar 

Method 

in  Com- 

parison  with  Experimental  and  Analytical  Results  Recorded  in  the 

Literature 


g 

P,IP^ 

_ 

5 

Pi  ' 

Literature 

i  Vector 

Literature 

Vector 

oc 

Expr. 

1  Analyl. 

'  Polar 

F.xper 

Xoalyt. ' 

Polar 

(/) 

3.0 

1  96 

1.9 

1  84 

2.3 

2.4 

2.18 

c 

u 

6.0 

2.9 

29 

2.86 

46 

49 

4.9 

ki. 

9.0 

3  45 

3  6 

3  63 

6  5 

1 

7.4 

1 

7.5 

£ 

C* 

5! 

0 

0 

2.15  (Ml  1.58) 

1.28 

— 

1.258 

261 

— 

2  66 

6.0(4/,  -  2.3) 

2.18 

^2  095 

5  43 

5.86 

V 

— 

X 

8 

— 

1  14.05 

14  2 

•a 

1 

c 

u 

10 

1  _ 

20.15 

20.0 

o 

X 

C. 

20 

— 

— 

70 

1  70 
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